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It  was  concluded  from  the  results  (an  ideal  control  system  had  been 
assumed)  that  it  would  be  most  difficult  to  design  a  pitch  control 
system  for  the  elimination  of  all  harmonic  airloads,  because  the 
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study  of  the  elimination  of  all  harmonic  airloads  and  r f  all  harmonic 
root  shears  is  not  planned. 

The  required  pitch  angle  inputs  applied  solely  at  the  blade  root  were 
computed  to  eliminate  the  noncancelling  oscillatory  root  shears.  Re¬ 
sults  indicate  that,  for  the  UH-1  configuration  considered,  most  of 
the  noncancelling  harmonic  root  shears  could  be  eliminated  with  second 
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formance  . 
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SUMMARY 


A  study  is  made  of  the  possibility  of  using  higher  harmonic 
pitch  angle  inputs  to  eliminate  the  transmission  of  oscillatory  vertical 
forces  from  a  helicopter  rotor  to  its  driving  shaft.  The  aerodynamic 
loads  that  are  the  source  of  these  objectionable  forces  are  computed 
by  the  use  of  a  realistic  model  which  represents  the  rotor  blades  by 
bound  vorticity  distributions  and  the  wake  by  a  mesh  of  segmented 
vortex  filaments. 

The  required  pitch  angles  are  first  found  which  will  eliminate 
the  oscillatory  lift  loadings  at  all  radial  stations  by  compensating  for 
the  azimuthal  variation  in  the  blade  tangential  velocity  and  for  the 
nonuniformity  of  the  downwash  in  the  rotor  disk.  In  the  computations 
for  this  case,  an  ideal  control  system  is  assumed  such  that  continuous 
radial  and  azimuthal  variations  in  pitch  angle  inputs  are  possible. 

In  the  second  part  of  the  study,  a  method  is  developed  for 
computing  the  pitch  angle  inputs  at  the  blade  root  required  to  eliminate 
the  oscillatory  root  shears  at  the  blade  attachment  hinge  without  making 
the  oscillatory  lift  equal  to  zero  at  all  radial  stations.  The  inertia  forces 
associated  with  blade  dynamic  motions  produce  root  shears  which  are 
taken  into  account  in  the  computations.  The  computational  procedure 
which  is  developed  is  capable  of  treating  either  the  two-bladed  teetering 
rotor  or  the  three-bladed  articulated  rotor. 

The  computed  results  are  based  on  a  two-bladed  teetering  rotor 
which  is  approximately  the  same  as  that  of  the  UH-1A  configuration  ex¬ 
cept  for  the  assumed  differences  in  pitch  control.  The  required  pitch 
angle  inputs  are  determined  for  several  cases  in  which  different  com¬ 
binations  of  harmonic  root  shears  are  eliminated. 

It  is  found  that  blade  dynamic  motions  have  an  important  effect 
on  the  pitch  inputs  required.  The  results  indicate  that  for  the  config¬ 
uration  considered,  most  of  the  oscillatory  vertical  root  shear  could 
be  eliminated  with  blade  root  pitch  inputs  at  the  second  and  fourth 
harmonics.  Considerable  variations  in  the  required  pitch  inputs  are 
found  for  different  flight  conditions. 
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INTRODUCTION 


High  fuselage  vibration  levels  are  still  today  a  source  of 
passenger  discomfort  in  many  helicopters  and  in  some  cases  are 
of  sufficient  magnitude  to  interfere  with  the  performance  of  desired 
helicopter  missions.  It  is  generally  believed  that  the  primary  sources 
of  the  objectionable  vibrations  are  the  rotor  generated  forces  which 
are  transmitted  directly  to  the  shaft  and  fuselage  via  the  blade  root 
fittings  and  rotor  hub.  This  report  is  limited  to  considering  the 
reduction  of  forces  of  this  type  --  more  specifically,  the  vertically 
transmitted  vibratory  root  shear  forces.  It  is  recognized,  however, 
that  it  is  also  necessary  to  minimize  vibration  from  other  sources 
(e.  g.  ,  those  caused  by  the  engine,  transmission,  tail  rotor,  and  the 
aerodynamic  excitation  of  the  fuselage  by  the  rotor  downwash). 

Several  different  methods  of  attack  have  been  made  on  the 
problem  of  reducing  the  vibratory  shears  which  are  transmitted 
directly  from  the  rotor  to  the  fuselage.  Passive  systems  for  isolating 
the  rotor  head  or  rotor-engine-transmission  system  from  in-plane 
vibratory  shears  have  been  widely  used  for  many  years  (e.  g,  ,  those 
discussed  in  Reference  1).  Some  success  has  been  achieved  with  this 
approach  in  spite  of  the  difficulty  of  isolating  against  the  large -ampli¬ 
tude,  low-frequency  vibrations  which  are  characteristic  of  the  helicopter 
rotor  system.  Currently,  work  is  being  done  on  various  active  vibration 
isolators  and  force  balancing  systems  (References  2  and  3).  However, 
it  will  be  some  time  before  such  systems  can  be  evaluated  in  service 
on  several  different  configurations. 

Another  philosophy  of  vibration  reduction  is  based  on  the  idea 
of  eliminating  the  exciting  forces  in  the  rotor  system  rather  than  iso¬ 
lating  them  from  the  rest  of  the  helicopter.  This  concept  is  well  ac¬ 
cepted  as  regards  the  elimination  of  mechanical  unbalance  or  forces 
caused  by  dissymmetry  between  blades.  Attempts  have  also  been  made 
to  reduce  the  amplification  of  the  response  in  the  natural  bending  modes 
ox  blades  by  the  addition  of  appropriate  ballast  weight  to  shift  the  natural 
bending  frequencies  and/or  increase  the  effective  aerodynamic  damping 
(Reference  4).  Reductions  in  bending  vibrations  obtained  in  this  manner 
would,  at  the  same  time,  reduce  the  associated  root  shears  transmitted 
to  the  shaft.  Another  approach  has  involved  the  application  of  second 
harmonic  pitch  control  both  to  reduce  vibration  exciting  forces  and  to 
obtain  certain  performance  benefits  (Reference  5). 

In  this  report,  further  consideration  is  given  to  the  possibility 
of  reducing  vertical  rotor  shaking  forces  by  the  application  of  suitable 
blade  pitch  angle  control.  The  work  reported,  which  has  been  largely 
exploratory,  started  with  finding  the  requirements  for  eliminating  all 
aerodynamic  exciting  forces  with  an  assumed  ideal  control  system 
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such  that  continuous  radial  and  azimuthal  variations  in  pitch  angle 
inputs  would  be  possible.  However,  emphasis  has  been  placed  on 
studying  some  of  the  inherent  dynamic  problems  in  eliminating  root 
shears  with  pitch  inputs  applied  solely  at  the  blade  root.  The  study 
has  been  limited  to  the  application  of  pitch  control  for  the  reduction 
of  vibrations,  and  the  effect  of  such  pitch  changes  on  aerodynamic 
performance  remains  to  be  determined.  Such  an  evaluation  should, 
of  course,  include  an  estimate  of  the  weight  and  power  require¬ 
ments  of  the  controlling  mechanism. 

In  this  first  study  of  the  problem,  it  was  decided  to  consider 
only  lift  or  vertical  shear  loads  transmitted  to  the  fuselage.  It  is  well 
known  that  vibratory  drag  loads  generated  in  the  rotor  system  are  also 
responsible  for  serious  fuselage  vibrations,  and  further  work  is  re¬ 
quired  to  determine  if  they  too  can  be  reduced  by  blade  pitch  control. 

The  vibratory  aerodynamic  forces  acting  on  a  rotor  blade  are 
dependent  upon  the  induced  velocities  produced  bv  the  vorticity  in  the 
wake  of  the  rotor.  Careful  consideration  of  the  effect  of  the  wake  vor¬ 
ticity  is  necessary  if  a  realistic  estimate  is  to  be  made  of  the  blade 
pitch  angles  required  to  eliminate  vibratory  blade  loads  and  root  shears. 
A  summary  is  given  in  Section  1  of  the  method  of  treatment  which  has 
been  used  for  the  aerodynamics  of  the  rotor-wake  system. 

The  objective  of  the  investigation  was  to  determine  the  spanwise 
and  azimuthal  distribution  of  blade  geometric  angles  of  attack  required 
to  eliminate: 


A.  All  harmonic  air  loads. 

B.  All  transmissible  {noncanceling)  harmonic  root  shears. 

C.  All  harmonic  root  shears. 

The  requirements  for  eliminating  all  harmonic  air  loads  (Case  A) 
are  taken  up  in  Section  2.  The  purpose  of  the  study  discussed  in  this 
section  is  to  find  the  radial  and  azimuthal  variations  in  blade  pitch 
angles  required  to  obtain  a  blade  lift  loading  which  does  not  vary  with 
azimuth  position  and  thus  is  optimum  from  the  standpoint  of  avoiding 
oscillatory  root  shears  and  blade  vibrations. 

However,  the  elimination  of  all  noncanceling  harmonic  root 
shears  (Case  B)  can  be  accomplished  without  imposing  the  severe 
requirement  of  an  azimuthly  independent  blade  load  distribution.  It 
is  only  necessary  chat  the  integral  of  the  harmonic  lift  and  inertia  force 
loadings  over  the  entire  rotor  be  zero;  the  blade  pitch  input  require¬ 
ments  are  made  easier  for  this  case  by  the  fact  that  the  root  shears  of 
the  separate  blades  cancel  at  some  harmonics.  In  the  event  of  dissym¬ 
metry  between  blades  preventing  perfect  cancellations,  all  transmitted 
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loads  could  be  eliminated  by  requiring  that  all  harmonic  root  shears 
be  zero  for  each  blade  individually  (Case  C).  It  would  be  expected 
that  consideration  of  Case  C  would  not  be  necessary  for  a  helicopter 
rotor  built  with  reasonable  manufacturing  tolerances. 

The  radial  variations  of  pitch  angles  required  to  make  the 
harmonic  root  shears  zero  are  not  unique.  A  computational  procedure 
is  formulated  in  Section  3  for  treating  Cases  B  and  C  which  assumes 
that  pitch  angle  inputs  are  introduced  by  one  or  two  control  modes 
at  each  harmonic. 

A  successive  approximation  solution  is  indicated  in  Section  3 
for  obtaining  blade  loads,  blade  motions,  and  required  pitch  control 
inputs.  The  current  approximations  for  the  blade  air  loads,  computed 
by  the  procedux-e  of  Section  1,  are  used  as  inputs  for  dynamic  response 
computations  which  give  new  blade  motions  to  be  used  in  the  compu¬ 
tations  for  the  next  approximation  for  the  air  loads.  The  procedure 
provides  for  including  flapwise  and  torsional  blade  dynamics  but 
neglects  in-plane  dynamic  motions. 

In  order  to  aid  convergence  of  the  overall  solution,  the  blade 
air  loads, which  are  computed  from  the  current  approximation  for  the 
blade  motions,  are  not  used  directly  in  the  response  computations  for 
the  next  approximation  but  are,  in  effect,  combined  with  incremental 
air  loads  based  on  estimates  of  the  change  in  air  loads  from  the  current 
to  the  next  approximation.  The  increments  in  air  loads  from  one  iteration 
to  the  next  become  smaller  and  smaller  as  the  solution  converges,  so  that 
the  final  solution  does  not  depend  on  the  method  used  for  estimating  the 
incremental  air  loads. 

It  was  hoped  that  satisfactory  convergence  of  the  solution  could 
be  obtained  using  estimated  incremental  air  loads  based  on  quasi-static 
aerodynamic  coefficients  which  could  be  precomputed  before  the  solution 
of  the  main  problem.  However,  it  was  found  that  the  incremental  air 
loads  based  on  quasi-static  theory  were  not  sufficiently  accurate  for 
all  cases. 

In  general,  convergence  difficulties  were  encountered  when 
the  torsional  degrees  of  freedom  were  included  in  the  computations. 

For  example,  a  torsional  mode  near  resonance  caused  convergence 
difficulties  in  computations  carried  out  for  the  reference  case  of 
conventional  helicopter  control  using  zero-harmonic  collective  and 
first-harmonic  cyclic  pitch  inputs.  The  computer  solution  for  this 
case  did  not  converge  when  the  estimated  incremental  air  loads  for 
this  torsion  mode  were  based  on  quasi-static  aerodynamic  theory. 
However,  satisfactory  convergence  was  obtained  when  an  improved 
estimate  of  the  incremental  torsional  air  loads  was  used  which  was 
based  on  coefficients  derived  from  the  change  in  loads  from  one 
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iteration  to  the  next  in  the  nonconvergent  solution.  These  same 
coefficients  were  tried  in  a  solution  for  Case  B,but  convergence 
difficulties  were  again  encountered  when  the  torsional  degrees  of 
freedom  were  included.  It  appeared  chat  further  refinement  of  the 
prediction  of  the  incremental  torsional  air  loads  would  be  necessary 
in  treating  Cases  B  and  C,but  this  could  not  be  attempted  within  the 
limits  of  the  research  program. 

Solutions  for  Case  B  which  were  obtained  by  a&suming  that  the 
torsional  deflections  were  zero  are  discussed  in  Section  4.  The  require¬ 
ments  for  eliminating  only  the  second-  and  fourth-harmonic  noncan¬ 
celing  shears  are  considered  as  well  as  those  for  the  elimination  of 
all  noncanceling  root  shears,  A  comparison  of  these  results  indicates 
some  of  the  problems  introduced  by  blade  flapwise  dynamics  in  elim¬ 
inating  oscillatory  root  shears  by  pitch  control. 

No  computational  results  are  presented  for  the  pitch  inputs 
required  to  eliminate  all  harmonic  root  shears  for  each  blade  individually 
(Case  C).  Convergence  difficulties  were  encountered  for  this  case  even 
when  the  torsional  deflections  were  set  equal  to  zero.  Methods  for  im¬ 
proving  the  estimated  incremental  air  load  in  order  to  obtain  conver¬ 
gence  were  not  studied  for  this  case. 

In  making  numerical  computations,  a  specific  configuration  has 
been  assumed.  Also,  several  assumptions  have  been  made  to  simplify 
the  computations  which  were  not  essential  to  the  methods  developed. 
Furthermore,  criteria  other  than  the  elimination  of  oscillatory  vertical 
lift  load  and  root  shears  have  not  been  used  in  finding  the  required  higher 
harmonic  blade  pitch  inputs.  In  order  to  clarify  the  basis  for  the  nu¬ 
merical  results,  the  principal  assumptions  and  limitations  made  in 
obtaining  them  are  summarized  below: 

1.  The  blade  feathering  and  torsional  elastic  axes 
are  assumed  to  be  coincident  at  the  quarterchord 
and  to  intersect  with  the  shaft  axis. 

2.  A  preconing  angle  of  3.  0  degrees  is  assumed. 

3.  The  blade  is  assumed  to  be  mass  balanced  about 
the  elastic  axis. 

4.  Compressibility  effects  are  not  considered. 

5.  The  wake  vorticity  is  represented  by  a  mesh  of 
vortex  filaments. 

6.  The  chordwise  vorticity  distribution  at  a  given 
radial  station  is  represented  by  a  Glauert  series. 
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7.  Three-dimensional  aerodynamic  effects  are  intro¬ 
duced  by  satisfying  normal  flow  boundary  conditions 
at  each  radial  station  including  the  induced  velocities 
due  to  wake  vorticity. 

8.  In-plane  root  shears,  flexibility,  and  deflections  are 
not  considered. 

9.  Rotor  performance  and  aircraft  flying  qualities  are 
not  considered  in  finding  required  pitch  angles. 

Additional  special  assumptions  which  are  applicable  to 
Sections  2  and  3  are  listed  at  the  beginning  of  those  sections. 
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1. 


AERODYNAMICS  OF  ROTOR-WAKE  SYSTEM 


The  treatment  of  the  aerodynamics  of  the  rotor-wake  system 
is  basically  the  same  as  that  used  in  Reference  6.  Consequently,  only 
an  abbreviated  description  of  the  aerodynamic  analysis  is  given  here, 
and  the  reader  is  referred  to  the  cited  reference  for  further  details 
and  a  justification  of  the  wake  model  used. 

The  shed  and  trailing  vorticity  distributions  in  the  wake  of 
each  blade  are  represented  by  an  arrangement  of  straight-line  vortex 
filaments  as  indicated  in  Figure  1.  For  simplicity,  the  case  of  only 


vf 


Figure  1  EXAMPLE  OF  WAKE  CONFIGURATION 
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four  radial  blade  segments  is  shown  on  the  diagram  i.  e.  ,  NR  =  4). 

Trailing  vortex  filaments  arise  only  at  the  ends  of  tie  blade  segments, 
since  it  is  assumed  that  there  is  no  radial  variation  in  bound  vorticity 
along  each  segment.  Shed  vorticity  is  deposited  continuously  in  the 
wake  behind  each  blade  of  an  actual  rotor  with  an  intensity  equal  to 
the  rate  of  change  of  bound  vorticity.  However,  the  numerical  solution 
is  carried  out  discontinuously  using  discrete  time  or  azimuth  angle 
increments.  The  discontinuous  azimuth  angles  are  denoted  by  V  with 
a  small  script  letter  used  as  a  subscript  (e.  g. ,  <p ^  or  ^  ).  The  equally 
spaced  azimuth  angles  are  given  by  the  formula  ^  -  2?r(A-/)/(NA) , 
where  A  takes  on  all  integer  values  from  one  to  the  total  number  of 
azimuth  positions  (NA).  The  shed  vorticity  deposited  by  each  blade 
segment  in  a  given  time  increment  is  lumped  into  a  single  shed  vortex 
filament  in  the  computational  model  as  indicated  on  the  sketch.  This 
concentrated  vortex  filament  is  shed  from  the  trailing  edge  at  the  posi¬ 
tion  corresponding  to  70  percent  of  the  time  increment.  Results  in 
Reference  6  show  that  this  procedure  gives  a  reasonable  approximation 
of  the  shed  wake. 

The  grid  of  straight-line  vortex  filaments  representing  the 
shed  and  trailing  vortex  filaments  can  be  truncated  after  a  prescribed 
number  of  azimuth  increments.  Further  aft,  the  wake  is  continued 
as  a  tip  trailing  vortex  filament  representing  the  rolled-up  vortex 
sheet.  Although  the  computational  program  permits  the  use  of  a 
distorted  wake,  the  motions  of  the  end  points  of  each  wake  segment 
were  computed  using  uniform  inflow  theory. 

In  carrying  out  the  solution,  boundary  conditions  are  satisfied 
for  NA  equally  spaced  azimuth  positions  of  the  rotor  at  the  midpoints 
of  the  NR  spanwise  segments  of  a  given  blade.  The  total  bound  vor¬ 
ticity  or  circulation  around  the  airfoil  for  each  of  these  (NR  x  NA  =•  NR  A) 
blade  segment  collocation  positions  is  denoted  by  r  with  a  lower  case 
letter  as  a  subscript  (e.  g.  ,  Rj  or  ).  The  integer  subscripts  j  or  t 
are  used  to  specify  successive  blade  segment  collocation  positions  for 
the  entire  rotor  disk  running  from  the  inboard  to  the  outboard  segments 
at  the  aft  azimuth  position  and  then  to  the  inboard  segment  at  the  next 
azimuth  position,  etc. 

For  a  periodic  problem,  the  strength  of  each  trailing  or  shed 
vortex  filament  can  be  expressed  by  a  linear  combination  of  the  fjs 
in  the  rotor  disk  as  a  consequence  of  the  vorticity  conservation  laws. 

When  the  bound  vorticity  strengths  in  adjacent  blade  segments  are  P 
and  r  +  t  ,  the  strength  of  the  trailing  vortex  segment  arising  at  their 
intersection  must  be  r  -PJ  +  f.  The  strength  of  the  shed  vortex  fila¬ 
ment  immediately  aft  of  a  given  blade  segment  is  equal  to  the  difference 
of  the  bound  vorticities  of  the  blade  segment  at  the  preceding  and  current 
azimuth  positions.  Similar  considerations  can  be  used  to  find  the  strengths 
of  all  the  vortex  filaments  in  a  wake  representation  such  as  shown  in 
Figure  1. 
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The  strength  of  the  trailing  vortex  segment  produced  by  the 
rolling  up  of  the  vorticity  deposited  in  the  wake  between  azimuth  posi¬ 
tions  and  is  denoted  by  ^  .  It  is  assumed  that  is  equal  in 

magnitude  to  the  largest  of  the  bound  vorticities  (T/s)  at  azimuth 
position  <p-k  •  This  amounts  to  assuming  that  all  the  trailing  vorticity 
deposited  in  the  wake  outboard  of  the  radial  position  of  maximum  bound 
vorticity  is  concentrated  in  the  tip  vortex. 

Each  segment  of  the  blade  at  each  discrete  azimuth  position 
is  represented  by  a  continuous  distribution  of  vorticity  aligned  in  the 
radial  direction  and  varying  in  intensity  only  in  the  chordwise  direction. 
Usirg  the  notation  of  Reference  6  as  indicated  in  Figure  2,  the  bound 
vorticity  of  a  blade  segment  at  collocation  position  (k)  can  be  expressed 
by  the  Glauert  series, 

[•  ^  ©O  _ 

A0k  coi-j  An^stnn9^  (1) 

n  z  1 

where 

X  -  -  b  cos  & .  (2) 


+  v 


Figure  2  REPRESENTATION  OF  BLADE  SECTION  BY 

CHORDWISE  BOUND  VORTICITY  DISTRIBUTION 
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The  condition  for  no  normal  flow  through  the  blade  chord  at 
the  midspan  of  a  particular  blade  segment  can  be  expressed  as  follows: 


( X )  *-  ( x) 


J_  r>b  yk  (f)d£ 
ZrrJ.b  *-<f 


=  o, 


(3) 


where 


and 


v'k  (»)  ~  (-h  +  V,os,9  +  6c9  *)k 

MRA 

u/,k(*) =  £  ck  (x)  rj . 

j  *  /  j 


(4) 

(5) 


The  symbol  v'k  denotes  the  velocity  relative  to  the  surface  of  the  airfoil 
due  to  plunging  motion  (Ti^Jand  geometric  angle  of  attack  f<%9k) while  u ^ 
is  the  normal  velocity  induced  by  the  trailing  and  shed  vortex  filaments 
in  the  wake.  It  is  noted  that  the  expression  for  «/*/.  can  be  written  in 
terms  of  the  bound  vorticity  strengths  because  the  strength  of  each 
wake  vortex  segment  is  a  function  of  the  ^  S  .  The  integral  in  Equa¬ 
tion  (3)  gives  the  induced  velocity  due  to  the  bound  vorticity.  An 
approximation  for  this  induced  velocity  is  used  which  should  be 
reasonable  at  the  midspan  of  the  k**  blade  segment;  it  is  based  on 
the  assumption  that  the  bound  vorticity  representing  the  blade 
segment  is  extended  to  infinity  in  both  directions.  This  explains  the 
use  of  the  tv/o -dimensional  Biot-Savart  law  in  the  integral  term. 

It  is  convenient  to  replace  the  chordwise  variable  x  in  Equa¬ 
tions  (3),  (4),  and  (5)  by  the  expression  defined  in  Equation  (2).  Then, 
the  evaluation  of  the  integral  in  Equation  (3)  gives  the  well-known  result, 

i  ft*  7ji(4)d<S  _/  r7r  yk(4>) sen  <pd4> 

ZrrJ^  x-  £  (c.os0  -  cos  0) 


oo 

=  Aok  "X  a»l  ^snQt  (6) 

n=  1 

where  the  series  in  Equation  (1)  has  been  substituted  for  Yk(<f>)  replacing 
0  by  j>  .  It  is  also  possible  to  expand  in  a  cosine  series  which  is 
written  in  the  form, 


bJRA  cq  oo 

(t»  *Z  (4.  *  2_  V  +Z  (r  ;Z  T  cos  ng)Q.  (7) 

J*t  '  J  n=>  J  /  V--/V  kl  nxl  kl  / 
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The  following  set  of  equations  results  from  substituting  Equa¬ 
tions  (4),  (6),  and  (7)  into  Equation  (3),  and  from  requiring  that  the 
coefficients  of  each  harmonic  cosine  term  be  zero  separately: 


A  relationship  between  the  r's  and  As  is  required  before  a  solution 
can  be  found.  This  is  obtained  by  integrating  the  Glauert  series  for 
the  bound  vorticity  over  the  chord  giving  the  following  expression  for 
the  total  bound  vorticity  or  circulation  about  blade  segment  k, , 

r” 

r ^  =  b  y^(9)  sin  OcLQ 


Since  the  theoretical  circulation  and  the  theoretical  lift  curve  slope 
are  not  achieved  in  practice,  an  empirical  correction  of  Equation  (9) 
is  made  by  multiplying  by  c^^jZrr,  giving 

b<‘(A»*TA')k-  (1° 

The  result  of  combining  Equations  (10)  and  (8)  is 
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where 


Tt  ■  (%)k  ^  (-  +  J  a, ,)k  , 

^.i '  fadt  ^  (\  "  i  -\ )  ’ 

and 

=  (C-*cc)k  bk  (T°kU  "  1r,kj)  * 


(12) 

(13) 


(14) 


Equation  (11)  represents  a  set  of  ( /V^/j^equations  which  are  solved  iter¬ 
atively  for  the  A/RA  fjs  by  the  Gauss -Sidel  method.  In  each  iteration, 
the  fj’s  for  the  different  spanwise  blade  segments  at  a  given  azimuth 
angle  (<fo)  must  be  compared  to  find  their  maximum  value  . 

The  tredtment  of  stalled  and  reversed  flow  regions  which  is 
used  in  Sections  3  and  4  is  carried  out  in  the  manner  described  in 
Reference  6.  It  is  assumed  that  stall  limits  the  maximum  value  of 
r-  attainable  at  each  blade  segment  collocation  position,  and  this 
maximum  value  is  introduced  as  an  additional  constraint  in  the  iter¬ 
ative  solution. 

Once  the  set  of  equations  for  the  /J's  (i.  e.  ,  Equation  (11))  has 
been  solved,  the  strengths  of  the  bound  vortices  of  the  blade  segments 
are  known  for  all  collocation  positions  in  the  rotor  disk.  Equation  (8) 
can  be  used  to  compute  the  corresponding  coefficients  (A^'s)  of  the 
Glauert  expansions  of  the  chordwise  vorticity  distributions  at  these 
positions.  The  time  derivatives  of  the  Glauert  coefficients  for  a  given 
blade  segment  are  also  required  in  the  computations.  They  are  deter¬ 
mined  by  assuming  that  all  variables  change  periodically  at  steady  flight 
conditions.  This  assumption  makes  it  possible  to  express  the  time 
derivatives  of  the  Glauert  coefficients  for  a  given  blade  segment  at  a 
particular  azimuth  position  in  terms  of  the  values  of  the  Glauert  co¬ 
efficients  for  the  blade  segment  at  all  of  the  Aft  equally  spaced  azimuth 
positions  which  are  used  in  the  computations. 

The  linearized  Bernoulli  equation  for  unsteady  flow  leads  to 
an  expression  for  the  chordwise  variation  in  pressure  difference  on  a 
blade  segment  as  it  moves  through  collocation  position  (b)  in  terms  of 
the  Glauert  coefficients  and  time  rate  of  change  of  the  Glauert  coefficients 
for  the  blade  segment  at  position  (k)  .  This  pressure  distribution  is  used 
to  compute  the  lift  and  pitching  moment  per  unit  span  as  a  function  of  the 
values  and  time  derivatives  of  the  Glauert  coefficients  of  a  blade  segment 
at  position  (k).  An  empirical  correction  factor  is  also  introduced  in 
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these  expressions  to  allow  for  the  fact  that  the  actual  circulation  is 
less  than  the  theoretical  value.  The  expressions  for  lift  and  pitching 
moment  (about  the  midchord)  per  unit  span  of  a  blade  segment  at 
collocation  position  ( k )  can  then  be  written  in  the  form, 

A  =  {cJctt’/,Yi(A0+{Al)  +  irb/>jj-(3Ac+A,  (15) 


-  {  J  (-/l,  +Az)-j  TTbf-  (A0+~At-  .(15) 

The  lift  and  moment  loadings  obtained  from  Equations  (15)  and 
(16)  are  those  re  alting  from  the  given  geometric  angles  of  attack 
rates  of  change  of  geometric  angle  of  attack (**9k)  ,  and  plunging  veloc¬ 
ities  (&k)  entering  in  Equations  (8)  and  (12).  These  quantities  must  be 
defined  in  terms  of  the  geometry  and  motion  of  the  rotor  system. 

The  shaft -oriented  reference  system  used  in  the  analysis  is 
shown  in  Figure  3  and  the  deflected  blade  axis  is  indicated  schematically 
in  Figure  4.  The  following  expressions  for  the  velocity  components  used 
in  the  analysis  are  readily  derived  by  reference  to  these  figures. 


Vf  =  fir  +  ( Vf  cos  cx,^)  -sin  ifj  (17) 

=  component  of  relative  velocity  perpendicular 
to  the  shaft  and  to  the  blade  axis. 

h  =  Vf  sln*scos(^j-bVf  cos  as  COS  Ip  sin  COS(^r) 

-  component  of  velocity  perpendicular  to  Vt  and 
the  blade  axis  (i.  e. ,  velocity  of  airfoil  relative 
to  air  excluding  induced  velocity). 

By  making  use  of  small  angle  approximations, 

.  JL 

h  -  Vf  sin  ccs  +  Vf  cos  cxs  cos  — 01  +  hm.  (18) 

dr 


Forward  velocity  is  denoted  by  V £,  shaft  angle  by  as^and  the  distance 
of  the  midchord  axis  above  the  reference  plane  by  h  . 
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Figure  3  RELATIVE  VELOCITIES  DUE  TO  FORWARD  MOTION  AND  BLADE 
ROTATION  (REFERRED  TO  A  SHAFT-ORIENTED  REFERENCE  SYSTEM) 


Figure  4  TRANSVERSE  VELOCITY  OF  BLADE  MIDCHORD  AXIS 
RELATIVE  TO  AIR  (INDUCED  YELOCIIlES  NOT  INCLUDED) 


2. 


REQUIREMENTS  FOR  ELIMINATING  ALL  HARMONIC  AIR  LOADS 


This  section  takes  up  the  requirements  on  the  spanwise  and 
azimuthal  variations  of  blade  geometric  angles  of  attack  which  are 
necessary  to  eliminate  all  harmonic  lift  loads  on  helicopter  rotors 
(Case  A),  and  thus  to  obtain  a  spanwise  lift  distribution  which  does 
not  vary  with  azimuth  position.  The  elimination  of  all  harmonic  loads 
will  be  defined  as  the  elimination  of  all  oscillatory  lift  loads  through 
the  twelfth  harmonic.  The  following  special  assumptions  are  made  in 
this  section  in  developing  a  method  for  finding  the  required  pitch  angles: 

1.  The  harmonic  lift  loads  at  all  spanwise  stations 
are  assumed  to  be  zero  for  the  first  through  the 
twelith  harmonic. 


2.  The  rotor  blade  is  assumed  to  be  modified  such  as  to 
preclude  the  reversed  flow  region. 

3.  Simplifications  are  made  in  the  Section  1  treatment 
of  rotor  aerodynamics  on  the  basis  of  assuming  the 
reduced  frequencies  to  be  small  at  all  harmonics 
of  importance. 


It  is  not  anticipated  that  the  pitch  change  schedule  which  is 
required  for  obtaining  the  desired  load  distribution  could  be  imple¬ 
mented  exactly  with  an  actual  helicopter.  However,  study  of  this 
ideal  case  is  useful  for  seeing  the  possibilities  and  limitations  of 
practical  systems. 

In  the  usual  method  for  the  computation  of  blade  loads,  including 
the  effect  of  wake  vorticity,  the  lift  is  determined  for  given  blade  pitch 
angle  changes.  This  procedure  might  be  iterated  several  times  with  the 
objective  of  finding  the  required  pitch  angles  to  eliminate  all  harmonic 
air  loads.  Instead,  the  inverse  problem  is  considered  which  consists 
of  solving  directly  for  the  geometric  blade  angles  necessary  to  produce 
a  constant  specified  lift  distribution.  Simplifying  assumptions  2  and  3, 
which  are  given  above,  are  made  in  order  to  make  practicable  the  solu¬ 
tion  of  the  inverse  problem.  These  approximations  are  believed  to  be 
sufficiently  accurate  for  a  general  study  of  the  problem. 


DISCUSSION  OF  SIMPLIFYING  ASSUMPTIONS  MADE  IN  ANALYSIS 


In  order  to  eliminate  nonlinearities  associated  with  reverse 
flow,  the  innermost  blade  airfoil  section  is  assumed  to  be  sufficiently 
far  outboard  to  eliminate  the  reverse  flow  region  completely;  any 
aerodynamic  forces  on  the  root  attachment  fitting  are  neglected. 


H 


A  comparison  is  made  in  Figure  5  of  the  standard  UH-1A  blade  plan- 
form  used  in  the  computations  of  Section  4  and  the  modified  blade  on 
which  the  computational  results  of  this  section  are  based.  The  blade 
feathering  and  elastic  axes  are  assumed  to  be  coincident  at  the  quarter  - 
chord  as  shown  in  Figure  5c. 

Since  the  aerodynamic  lift  distribution  for  Case  A  is  assumed 
to  be  constant  as  the  rotor  turns,  there  will  be  no  aerodynamic  forces 
tending  to  excite  oscillatory  teetering  or  bending  motions  of  the  blades 
relative  to  the  rotor  disk.  If,  in  addition,  the  blade  is  assumed  to  be 
dynamically  balanced,  there  is  no  need  for  considering  blade  bending 
dynamics  in  the  analysis. 

The  computations  for  the  required  spanwise  and  azimuthal 
variations  in  blade  geometric  angles  of  attack  give  the  sum  of  the  pitch 
control  and  torsional  deflections  at  each  radial  station.  Since  no  con¬ 
trol  system  is  being  postulated  here,  it  is  immaterial  how  the  total 
blade  angle  is  divided  between  pitching  and  torsional  motions.  Tor¬ 
sional  dynamics  are  not  included  in  the  analysis  because  calculations 
are  not  performed  to  separate  the  total  angles  into  torsional  deformations 
and  control  inputs.  Meaningful  computations  for  torsional  deformations 
cannot  be  made  without  specifying  a  control  input  mechanism. 

In  order  to  indicate  other  simplifying  assumptions  which  are 
made,  it  is  convenient  to  introduce  the  dimensionless  time, 

and  rewrite  Equations  (10),  (11),  and  (15)  in  the  form 

tjkh*(Ao  +  iAt)k  •  (i9) 


rk 


UbooA 

vi  j 


k 


V/?/?  HA 

+  H  °~k-rj  r^'  ^  •  (20) 

/  J  *  -  ■  At 

J  =  '  J  -  / 
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(USED  IN  COMPUTATIONS  OF  SECTION  4) 


c.  ASSUMED  LOCATION  OF  BLADE  PITCHING  AND  ELASTIC  AXES 


Figure  5  BLADE  SEGMENTS  USED  IN  COMPUTATIONS 
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and 


4  ■  ^  ^  {h(***iAbi  (7* )k'  (JA°+Ai+ 7A*)\k- 


The  magnitude  of  the  characteristic  frequency  u)c  is  to  be 
selected  such  that  the  derivatives  of  the  variables  with  respect  to 
dimensionless  time,  t  ,  are  of  the  same  order  or  smaller  than  the 
original  variables.  It  is  clear  that  the  characteristic  frequency  will 
be  of  the  order  of  the  highest  harmonics  of  importance  in  a  Fourier 
series  description  of  the  loading.  The  simplifying  assumption  is  made 
that 


ba>c 

—  «  /,  (22) 


which  is  equivalent  to  assuming  that  the  reduced  frequencies  at  all 
harmonics  of  importance  are  small.  On  the  basis  of  the  assumption 
in  Equation  (22),  the  term  in  Equation  (20)  and  the  time  deri¬ 

vative  term  in  Equation  (21)  are  dropped.  In  Equation  (20),  the  dimen¬ 
sionless  linear  deflection  (b/r)  is  of  the  order  of  the  variable  part  of  ccg 
in  the  outer  part  of  the  blade.  On  the  other  hand,  the  ratio  of  the  radial 
position  to  the  semichord,  (f/b)  ,  becomes  a  large  number  near  the  blade 
tip.  Thus,  it  appears  that  the (d/dt )(h/r) term  in  Equation  (20)  may  be  of 
larger  magnitude  than  the  other  time  derivative  term  and  is  retained. 
When  the  approximations  indicated  above  are  made  and  Equation  (19) 
is  substituted  into  Equation  (21),  the  following  simplified  equations 
result  which  are  the  basis  for  the  analysis  in  this  section: 

a'/?/?  ha 

rk  -  cj  bk  (V,  a9-h)k  +Z.  ^k:  rj  z*  r2  (23) 


•A  =  [pVf  r)k  . 


(24) 
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Equation  (18)  is  u»ed  to  evaluate  h  -  dhjdt  setting 


and 


dr 


•JO* 


dhQ 

dr 


•> 


/B  =  /Bfc  COS  tp  +  J0{s  sin  p, 


where  fi  is  the  first  harmonic  flapping  relative  to  a  plane  perpendicular 
to  the  shaft  and  ha  is  the  deflection  due  to  preconing  and  steady  blade 
bending.  The  term  in  hm  is  dropped  in  accordance  with  the  assump¬ 
tions  made  in  simplifying  Equation  (20).  After  substituting  the  above 
relations  into  Equation  (18),  the  following  expression  for  H  results: 

h  -  Vf  sin  a-s  +■  Vf  cos  ccs  cos  tp 

■t  (-  r£ifiic  sin  tp  +  rnpu  cos  p).  (25) 


(~f-  +fitc  cos  +Ais  ?) 


METHOD  OF  SOLUTION  OF  INVERSE  PROBLEM 


As  mentioned  previously,  a  constant  lift  distribution  along  the 
blade  is  specified  in  the  solution  of  the  "inverse  blade  loading  problem”. 
A  number  of  possible  criteria  might  be  considered  in  specifying  the 
loading  distribution  (e.  g.  ,  one  to  obtain  maximum  lift,  minimum  power, 
etc. ).  In  the  present  program,  it  was  desired  that  computations  be 
carried  out  for  typical  rotor  operating  conditions  for  which  data  were 
available.  The  conditions  selected  for  analysis  were  approximately 
the  same  as  those  obtained  in  flight  test  conditions  DN66,  DN67,  and 
DN68  of  Reference  7  which  were  at  advance  ratios  of  A-  =  0.  215,  0.  259 
and  0.251,  respectively. 

Figure  6  presents  the  distribution  of  effective  angle  of  attack 
of  the  standard  UH-1A  configuration  at  p  -  0.  259  as  obtained  by  compu¬ 
tations  discussed  in  Section  4.  The  reverse-flow  region  is  crosshatched 
while  the  boundary  of  the  stall  region  is  shown  with  a  dashed  line.  The 
assumptions  lor  eliminating  the  reverse-flow  region  have  already  been 
discussed.  Elimination  of  the  stall  region  can  be  accomplished  by  de¬ 
creasing  the  effective  angle  of  attack  in  the  inboard  portion  of  the  span 
but  compensating  increases  must  be  made  in  the  effective  angles  of 
attack  near  the  blade  tip  if  the  total  thrust  is  to  be  unchanged.  It  was 
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found  that,  in  order  to  obtain  the  thrust  corresponding  to  the  }J-  =  0.  259 
condition  with  an  azimuthally  independent  load  distribution,  the  retreating 
blade  would  have  to  be  operating  close  to  stall  at  all  radial  positions. 

Therefore,  it  is  assumed  that  the  specified  lift  distribution  is 
the  one  obtained  at  azimuth  angle  f  -  270  degrees  with  constant 
along  the  span; 

^  ~  cj  P  (^D)(Vt)  f-z7o  (26) 


=  J  P  (2b)(r±r  -  Vf  cos  asf.  (27) 


In  Figure  7,  the  lift  distribution  given  by  Equation  (27)  is  compared 
with  the  computed  lift  of  the  standard  UH-1A  blade  at  an  advance  ratio 
of  /4  =  0.  259.  The  magnitude  of  in  Equation  (27)  is  selected  such 
that  the  integrated  lift  loading  gives  the  desired  total  lift  per  blade. 

As  mentioned  previously,  the  Cj  required  for  the  jj.  -  0.  259  flight 
condition  results  in  the  blade's  being  close  to  stall  at  <f-'  =  270  degrees. 
More  freedom  in  rotor  design  would  be  desirable,  making  it  possible 
to  obtain  the  same  total  lift  while  operating  at  lower  section  lift  co¬ 
efficients.  This  might  be  accomplished  by  a  proper  compromise  of 
blade  taper  and  changes  in  average  chord  and  blade  radius.  If  the 
blade  design  were  revised  in  this  manner,  it  might  also  be  possible 
to  choose  a  specified  lift  distribution  with  a  lower  tip  loading  in  order 
to  reduce  compressibility  effects  on  the  advancing  blade.  Compressi¬ 
bility  effects  have  not  been  considered  in  this  work. 

The  lift  loading  specified  by  Equation  (27)  is  to  be  maintained 
constant  at  all  azimuth  positions  and  when  equated  to  Equation  (24) 
gives  the  following  expression  for  the  bound  vorticity  strength  at  an 
arbitrary  blade  segment  position  k  : 


r k.  =  clbk  (a  rk  -  Vf  cos  as)  cj V,k  .  (28) 
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Figure  7  COMPUTED  LIFT  LOADING  DISTRIBUTIONS  FOR  UH-1A  AT  y.  -  0.259 
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Having  found  the  bound  A s  for  all  radial  and  azimuth  positions  by 
Equation  (28),  the  geometric  angles  of  attack  required  to  obtain  the 
prescribed  loading  can  be  computed  from  Equation  (23)  rearranged 
as  follows: 


(X. 


3k 


_1_ _ 

to Jk  **  \ 


NA 

r‘ 


(29) 


* 

The  term  is  computed  from  Equation  (25)  for  a  prescribed  tip-path- 
plane  position  and  the  average  deflection  due  to  the  prescribed  loading. 


RESULTS  OF  COMPUTATIONS  FOR  PITCH  ANGLES  REQUIRED  TO 

Eliminate  all  haRmoNiC  Aik  LOADS 


Computations  were  carried  out  fo:.  three  different  flight  con¬ 
ditions  using  the  parameters  shown  in  Table  I,  These  conditions  are 
approximately  the  same  as  those  that  existed  in  flight  test  conditions 
DN66,  DN67,  and  DN68  of  Reference  7,  but  they  have  been  adjusted 
slightly  so  that  the  effects  of  varying  forward  speed  and  rotor  angular 
velocity  can  be  seen  independently.  The  adjusted  flight  conditions  are 
denoted  as  DN66A,  DN67A,  and  DN68A.  Table  I  presents  the  shaft 
angles  and  the  tip -path -plane  positions  required  for  trimmed  flight  at 
these  conditions  in  addition  to  the  advance  ratios,  forward  velocities, 
and  rotor  angular  velocities.  The  rotor  planform  used  in  the  compu¬ 
tations  differs  somewhat  from  the  UH-1A  configuration  because  of 
moving  the  blade  root  outboard,  as  discussed  previously. 

Figure  8  shows  azimuthal  variations  in  blade  pitch  angles 
(*'’tch  control  inputs  plus  torsional  deflections)  which  are  required  to 
e  .rninate  all  nonsteady  air  loads  (lift  loading)  through  the  twelfth  har¬ 
monic.  The  plots  which  are  shown  for  conditions  DN66A  and  DN67A 
include  only  the  variable  part  of  the  required  pitch  angles,  and  the  mean 
angles  required  are  given  for  each  radial  position.  The  forward  velocity 
for  condition  DN66A  is  17.  3  percent  less  than  the  forward  velocity  for 
condition  DN67A,but  the  rotor  angular  velocities  for  these  conditions 
are  the  same. 

Plots  of  the  variable  part  of  the  required  pitch  angles  are  not 
presented  for  flight  condition  DN68A.  The  angular  velocity  for  condition 
DN68A  is  3,  5  percent  greater  than  the  one  for  condition  DN67A,  while 
the  forward  velocities  for  the  two  cases  are  the  same.  The  azimuthal 
plots  for  condition  DN68A  did  not  depart  significantly  from  those  shown 
on  Figure  8  for  condition  DN67A  and  are  not  included  to  avoid  confusing 
the  curves.  However,  results  for  condition  DN68A,  which  are  included 
in  the  tabular  data  to  be  discussed,  indicate  differences  in  the  harmonic 
content  of  the  required  pitch  angles  for  conditions  DN67A  and  DN68A. 
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Table  1 

PARAMETERS  FOR  ADJUSTED  FLIGHT  CONDITIONS 


PARAMETERS 

ADJUSTED 

FLIGHT  CONDITION* 

DN66A 

DN67A 

DN68A 

V  =  FORWARD  VELOCITY  (ft/sec) 

(knots) 

155.4 

92. 1 

188.0 

1 11.4 

188.0 

III. 4 

IL  =  ROTOR  ANGULAR  VEL.  (rad/sec) 

( rpra) 

32.8 

313 

32.8 

313 

33.9 

324 

fX  =  ADVANCE  RATIO 

0.215 

0.259 

0.251 

/>  =  AIR  DENSITY  (slugs/ft3) 

0.002155 

0.002155 

0.002155 

K4  =  LOAD /BLADE  (lb) 

3200 

3200 

3200 

<XS  =  SHAFT  ANGLE  (deg) 

4.6 

6.5 

6.5 

TIP-PATH-PLANE  POSITION 

RELATIVE  TO  SHAFT: 

A  -  FORWARD  TILT  Idea) 

'-IC 

n  ee 

”  V  •  U  'J 

A  AP 

-v  .o J 

fi1s  =  LEFT  TILT  (deg) 

0.68 

1.80 

1.46 

‘the  ADJUSTED  FLIGHT  CONDITIONS  DN66A,  DN67A,  AND  DN68A  ARE  SLIGHTLY  DIFFERENT 
FROM  THE  CORRESPONDING  CONDITIONS  OF  REFERENCE  7. 
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PITCH  ANGLE  -  degrees  PITCH  ANGLE  -  degrees  PITCH  ANGLE  -  degrees  PITCH  ANGLE  -  degrees 


FLIGHT  CONDITION: 

-  0N66A 

Vf  =  155.1  ft/sec 
Cl  -  32.8  rad/sec 
M-  =  0.215 

-  DN67A 

Vf  =  188.0  ft/sec 
fl  =  32.8  rad/sec 
M-  =  0.259 


Pirtjiro  O 

•  ■  uw<  v  \J 


AZIMUTHAL  VARIATIONS  OF  BLADE  PITCH  ANGLES  REQUIRED  TO 
ELIMINATE  ALL  HARMONIC  AIR  LOADS  (MODIFIED  UH-1A  CONFIGURATION) 
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It  can  be  seen  from  the  plots  on  Figure  8  that  ordinary  first 
harmonic  pitching  motion  is  by  far  the  largest  harmonic  component 
present.  These  first  harmonic  pitch  changes  are  all  of  approximately 
the  same  magnitude  with  the  exception  of  those  at  the  most  inboard 
station;  they  result  primarily  from  the  requirements  on  rotor  tip-path- 
plane  position.  Second  and  higher  harmonic  components  can  be  observed 
in  the  azimuthal  variations  of  the  required  pitch  angles  but,  in  general, 
the  curves  appear  fairly  smooth  except  for  irregularities  found  at  the  most 
inboard  station. 

In  order  that  a  more  qr  antitative  evaluation  of  the  data  can  be 
made,  the  results  of  the  harmonic  analysis  of  the  required  blade  pitch 
angles  (i.  e. ,  the  required  geometric  angles  of  attack)  are  presente  I  in 
Table  II.  At  the  inboard  station,  r/R-  0.400,  where  large  forward  velocity 
variations  occur,  a  second  harmonic  pitch  angle  change  of  4.  2  degrees 
is  required  for  condition  DN67A.  The  required  second  harmonic  pitch 
variation  is  reduced  to  2.  1  degrees  at  the  next  outboard  station (r/R- 
0.  588).  For  further  outboard  stations,  the  required  second  harmonic 
control  inputs  do  not  exceed  1,  1  degrees. 

The  tables  indicate  that  the  required  pitch  control  inputs  at 
harmonics  from  the  third  to  the  fifth  are  of  the  order  of  0.  5  degree 
except  for  r/R-  0.  400  where  they  reach  1. 8  degrees.  The  required 
pitch  angles  for  the  sixth  to  the  twelfth  harmonics  are  considerably 
lower. 


Although  the  plots  for  flight  condition  DN67A  (  n  =  32.  8  rad/ sec) 
and  flight  condition  DN68A  ( -O.  =  33,  9  rad/sec)  were  almost  indistin¬ 
guishable,  it  can  be  seen  in  Table  II  that  small  differences  exist  in  both 
the  amplitudes  and  phases  of  the  pitch  angles  required  to  eliminate  the 
various  harmonic  air  loads.  Much  larger  differences  in  the  required 
amplitudes  and  phases  of  the  pitch  angles  are  found  for  the  speed  incre¬ 
ment  between  flight  condition  DN66A  (  Vf  =  155.  4  ft/sec)  and  flight  con¬ 
dition  DN67A  (  Vf.  -  188.0  ft/ sec). 

In  order  to  get  a  better  idea  of  how  the  ideal  control  angles 
computed  in  this  section  might  be  approximated  by  an  actual  control 
mechanism,  plots  are  presented  on  Figure  9  of  required  pitch  angles 
versus  radial  position.  The  curves  are  presented  for  increments  in 
azimuth  angle  of  30  degrees.  Similar  curves  were  also  obtained  for 
azimuth  positions  intermediate  to  the  ones  given  but  a^e  not  presented 
since  tney  show  the  same  trends  as  the  other  results.  In  Figure  9  as 
in  rigure  8,  results  are  given  corresponding  only  to  flight  conditions 
DN66A  and  DN67A,  since  the  results  for  DN68A  were  very  close  to 
those  for  DN67A. 
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Table  II 

HARMONIC  ANALYSIS  OF  GEOMETRIC  ANGLES  OF  ATTACK 
REQUIRED  TO  ELIMINATE  HARMONIC  AIR  LOADS 


(a)  r/R  =  0.400 


FLIGHT 

CONDITION 

DN66A 

DN67A 

DN68A 

HARMONIC 

AMPLITUDE 

PHASE 

AMPLITUDE 

PHASE 

AMPLITUDE 

PHASE 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

0 

G.856 

0.0 

8.842 

0.0 

6.348 

0.0 

1 

6.171 

40.6 

9.151 

54.1 

8.094 

53.7 

2 

2.573 

146.1 

4.208 

173.9 

3.724 

174.0 

3 

0.358 

33.8 

1.783 

-84.0 

1.484 

-79.5 

4 

0.916 

-178.1 

1.254 

-15.9 

0.988 

-24.0 

5 

1.192 

138.3 

0.648 

152.6 

0,702 

151.9 

6 

0.543 

-33.6 

0.455 

135.7 

0.361 

111.3 

7 

0.249 

-103.9 

0.172 

-42.8 

0.204 

-62.6 

8 

0.359 

32.6 

138.1 

0.136 

125.0 

9 

0.341 

-128.7 

42.1 

0.258 

18.2 

10 

0.486 

130.7 

0.175 

-60.8 

0.274 

-86.3 

0.213 

21.8 

0.123 

-147.7 

0.178 

165.8 

f 

0.102 

-180.0 

0.032 

-180.0 

0.086 

0.0 

(b)  r/R  =  0.588 


0 

8.652 

0.0 

9.548 

9.048 

0.0 

1 

5.856 

51.5 

7.419 

6.731 

57.4 

2 

1.056 

149.0 

2.095 

1.716 

-178.0 

3 

0.664 

114.1 

0.610 

0.365 

-99.3 

4 

0  243 

135.7 

0.367 

37.2 

55.8 

5 

0.579 

-118.7 

0.047 

-114.0 

0.0C8 

-179.2 

6 

0.148 

-28.6 

0.346 

-172.0 

0.359 

-178.7 

7 

0.662 

10.8 

0.271 

123.3 

0.246 

119.4 

8 

0.226 

-33.5 

0.246 

28.5 

0.229 

21.4 

9 

0.366 

172.2 

0.195 

-56.6 

0.193 

-57.3 

10 

0.113 

151.8 

0.168 

-143.4 

0.155 

-139.9 

11 

0.269 

-67.2 

0.104 

115.8 

0.082 

122.2 

12 

0.013 

0.0 

0.080 

0.0 

0.052 

0.0 

AZIMUTHAL  VARIATION  OF  REQUIRED  nth  HARMONIC  ANGLE  =  ( Amp  1 . )  x  cos  (nf  +  Phase  Angle) 
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Table  II 


HARMONIC  ANALYSIS  OF  GEOMETRIC  ANGLES  OF  ATTACK 
REQUIRED  TO  ELIMINATE  HARMONIC  AIR  LOADS 


(c)  r/R  =  0.750 


FLIGHT 

CONDITION 

DN66A 

DN67A 

0N68/ 

HARMONIC 

AMPLITUDE 

PHASE 

AMPLITUDE 

PHASE 

AMPLITUDE 

PHASE 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

0 

9.1111 

0.0 

10.630 

0.0 

9.995 

0.0 

1 

1.991 

58.1 

6.788 

59.0 

6.108 

58.6 

2 

0.366 

161.0 

1.130 

-178.1 

0.881 

177.9 

3 

0.170 

168.1 

0.273 

-71.6 

0.119 

-67.8 

1 

0. 408 

150.9 

0.381 

96.1 

0.392 

103.6 

5 

0.181 

-108.6 

0.259 

101.7 

0.223 

113.0 

6 

0.101 

-65.1 

0.082 

-151.3 

0.110 

-121.0 

7 

0.1711 

-10.1 

0.153 

-150.7 

0.131 

-126.2 

8 

0.198 

72.8 

0.083 

121.6 

0.061 

61.1 

9 

O.lli* 

83.5 

0.  i  67 

30.1 

0.205 

26.1 

10 

0.1 51 

-101.8 

0.177 

-21.8 

0.160 

-12.1 

0.232 

-112.7 

0.132 

-101.1 

0.116 

-118.5 

0.098 

-180.0 

0.151 

-180.0 

0.172 

-180.0 

(d)  r/R  =  0.850 

0 

10.100 

0.0 

11.390 

0.0 

10.660 

0.0 

1 

1.569 

60.1 

6.531 

59.5 

5.813 

59.1 

2 

0.300 

175.9 

0.810 

176.9 

0.661 

171.5 

3 

0.097 

-158.6 

0.226 

-56.7 

0.119 

-53.9 

1 

0.351 

153.5 

0.283 

108.8 

0.288 

1 18.5 

5 

0.132 

-117.7 

0.257 

121.9 

0.212 

133.5 

6 

0.301 

-71.6 

0.052 

159.1 

0.062 

-121.3 

7 

0.179 

-67.1 

0.095 

-88.7 

0.138 

-75.8 

8 

0.116 

71.6 

0.071 

-57.1 

0.101 

-36.1 

9 

0.175 

63.8 

0.075 

11.0 

0.115 

17.7 

10 

0.079 

-11.9 

0.071 

25.0 

0.079 

31.0 

II 

0.193 

-120.6 

0.009 

11.9 

0.015 

171.1 

12 

0.  198 

-180.0 

0.038 

-180.0 

0.113 

-180.0 

AZIMUTHAL  VARIATION  OF  REQUIRED  nth  HARMONIC  ANGLE  ^  (Ampi.)  x  cos  (nf  +  Phase  Angle) 
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Table  II 

HARMONIC  ANALYSIS  OF  GEOMETRIC  ANGLES  OF  ATTACK 
REQUIRED  TO  ELIMINATE  HARMONIC  AIR  LOADS 

(e)  r/R  =  0.900 


FLIGHT 

CONDITION 

DN66A 

DN67A 

DN68A 

HARMONIC 

AMPLITUDE 

PHASE 

AMPLITUDE 

PHASE 

AMPLITUDE 

PHASE 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

0 

0.0 

11.830 

0.0 

11.050 

■ESI 

1 

■mifl 

63.5 

6.396 

61.1 

5.671 

mum 

2 

BSKSb. 

-176.3 

0.733 

178.5 

0.571 

BBSS 

3 

0.153 

-63.0 

0.252 

-47.4 

0.202 

-46.0 

4 

0.347 

170.2 

0.231 

121.9 

0.232 

133.1 

5 

0.293 

177.1 

0.277 

135.1 

0.251 

144.7 

6 

0.104 

-84.3 

0.064 

133.8 

-166.5 

7 

0.240 

-38.9 

0.124 

-54.3 

0.165 

-52.5 

8 

0.158 

-28.4 

0.149 

-57.4 

0.171 

-48.0 

9 

0.044 

69.4 

0.050 

••31.4 

0.061 

-17.2 

>0 

0.110 

124.3 

0.076 

81.8 

0.081 

91.6 

11 

0.094 

138.6 

0.081 

86.8 

0.091 

117.7 

12 

0.066 

-180.0 

0.010 

0.0 

0.056 

-180.0 

0 

11.460 

0.0 

12.920 

0.0 

12.040 

0.0 

1 

4.496 

67.9 

6.484 

64.9 

5.734 

65.1 

2 

0. 189 

-14a. 0 

0,736 

- i 73 • i 

0,577 

;  -173,2 

3 

0.202 

-53.0 

0.268 

-42.2 

0.221 

-41.3 

4 

0.331 

170.1 

0.204 

12C.1 

0.207 

137.3 

5 

0.345 

168.2 

0.260 

140.0 

0.246 

149.0 

6 

0.052 

-135.6 

0.059 

141.9 

0.040 

-178.8 

7 

0.249 

-38.7 

0.135 

-52.0 

0.170 

-51.2 

8 

0.226 

-38.8 

0.167 

-58.2 

-52.1 

9 

0.021 

32.1 

0.051 

-60.8 

-38.4 

10 

0.183 

115.4 

0.076 

98.6 

0.091 

104.4 

11 

0,175 

113,4 

0.092 

93.0 

0.108 

112.8 

12 

0.047 

-180.0 

0.004 

0.0 

0.038 

-180.0 

AZIMUTHAL  VARIATION  OF  REQUIRED  nth  HARMONIC  ANGLE  =  (Ampl.)  x  cos  (ntp  +  Phase  Angle) 
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Table  II 

HARMONIC  ANALYSIS  OF  GEOMETRIC  ANGLES  OF  ATTACK 
REQUIRED  TO  ELIMINATE  HARMONIC  AIR  LOADS 

(g)  r/R  =  0.987 


FLIGHT 

CONDITION 

DN66A 

DN67A 

DN68A 

HARMONIC 

AMPLITUDE 

PHASE 

AMPLITUDE 

PHASE 

AMPLITUDE 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

0 

0.0 

14. 580 

0.0 

0.0 

1 

mw 

75.7 

6.905 

72.3 

72.7 

2 

-131.5 

0.920 

-159.5 

-158.6 

3 

-45.2 

0.337 

-36.3 

-36.0 

4 

166.0 

0.181 

122.1 

135.6 

5 

0.341 

166.6 

0.243 

145.0 

153.2 

6 

0.071 

-156.7 

0.068 

158.6 

0.058 

-170.0 

7 

0.222 

-41.3 

-55,9 

0.152 

-54.2 

8 

-37.0 

0.157 

-55.1 

0.178 

-49.8 

9 

19.8 

0.059 

-51.2 

0.061 

-39.0 

10 

109.1 

0.062 

97.6 

0.078 

104.4 

11 

116.3 

0.082 

102.5 

0.10! 

118.5 

12 

0.074 

-180.0 

0.018 

-180.0 

0.057 

-180.0 

AZIMUTHAL  VARIATION  OF  REQUIRED  nth  HARMONIC  ANGLE  =  (Ampl.)  x  cos  (n^  +  Phase  Angle) 


It  is  seen  from  the  curves  on  Figure  9  that  the  radial  variation 
of  required  pitch  angles  at  each  azimuth  position  can  be  fitted  quite 
well  with  one  or  two  straight-line  segments  in  the  case  of  the  higher 
speed  flight  condition  DN67A.  However,  such  a  representation  of  the 
curves  for  condition  DN66A  is  less  satisfactory  because  some  of  them 
exhibit  a  pronounced  "S"  shape. 

Figure  10  is  presented  in  order  to  show  the  mean  pitch  angle 
versus  radius  distributions  found  in  the  computations  for  the  required 
pitch  angles  to  eliminate  all  harmonic  air  loads.  It  can  be  seen  that 
the  appropriate  mean  twist  distribution  varies  with  the  flight  regime. 
Consequently,  designing  for  the  entire  flight  spectrum  would  necessitate 
selecting  a  built-in  blade  twist  distribution  which  would  be  a  compromise 
between  the  distributions  required  for  low-  and  high-speed  flight.  The 
increments  in  mean  twist  from  the  built-in  twist  as  required  for  different 
flight  conditions  would  place  additional  demands  on  the  blade  pitch  angle 
control  system. 


PITCH  ANGLE  -  degrees  PITCH  ANGLE  -  degrees  PITCH  ANGLE  -  degrees  PITCH  ANGLE  -  degrees 


<P  =  60° 


210 


f  -  240° 


0 

0.4  0.6  0.8  1.0 

O’  0.4  0.6  0.8 

1.0  0  0.4  0.6  0.8  1.0 

r/R 

r/R 

r/R 

FLIGHT  CONDITION: 

DN66A 

DN67A 

- - 

—  Vf  =  155.4  ft/sec 

-  Vf  =  188.0  ft/sec 

tl  -  32.8  rad/sec 

O  —  00  A  rorl  / 

=  0.215 

A-  =  0.259 

Figure 

9  RADIAL  VARIATION 

OF  BLADE  PITCH  ANGLES  REQUIRED  TO  ELIMINATE 

ALL  HARMONIC  AIR 

LOADS  (MODIFIED  UH-1A  CONFIGURATION) 
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Figure  10  REQUIRED  MEAN  PITCH  ANGLE  VERSUS  RADIUS  DISTRIBUTION 
TO  ELIMINATE  ALL  HARMONIC  AIR  LOADS 


The  large  positive  twist  distribution  indicated  on  Figure  10  is 
in  marked  contrast  with  present  design  practice  in  which  a  negative  twist 
distribution  is  used.  It  should  be  emphasized  that  the  required  positive 
twist  distributions  which  were  computed  are  a  direct  consequence  of  the 
specified  lift  distribution  shown  on  Figure  7.  This  lift  distribution  had 
an  unusually  high  tip  loading.  A  positive  mean  twist  distribution  is  not 
a  general  characteristic  of  the  pitch  angle  inputs  required  to  obtain  an 
azimuthally  independent  lift  distribution. 


f- 
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3.  METHOD  FOR  FINDjlNG  PITCH  ANGLE  INPUTS 

TO  ELI  MINA  TTTOSCILLATORY  ROOT  SHEAR'S* 


In  the  preceding  section,  the  blade  pitch  angle  variations  were 
found  which  are  necessary  to  eliminate  all  harmonic  lift  loads  on  a 
rotor  and,  consequently,  all  harmonic  shears  at  the  blade  root.  A 
method  is  considered  in  this  section  for  finding  the  pitch  angle  inputs 
to  satisfy  the  less  severe  requirements  for  eliminating  only  oscillatory 
vertical  root  shears.  The  continuous  radial  variation  of  pitch  angle 
control  which  was  used  in  Section  2  is  no  longer  assumed,  and  the  root 
shears  must  be  made  zero  without  making  all  harmonic  air  loads  zero. 
The  blade  dynamic  motions  excited  by  these  residual  harmonic  air  loads 
affect  the  root  shears  and  are  included  in  the  analysis. 

In  particular,  a  computational  procedure  is  developed  for  finding 
required  harmonic  control  inputs  in  the  case  of  the  two-blade  teetering 
rotor  assuming: 

1.  The  flight  condition  is  steady  so  that  blade  dynamic 
motions  can  be  assumed  periodic. 

2.  The  blade  root  shears  due  to  combined  lift  and 
inertia  force  loadings  are  to  be  zero  at  prescribed 
harmonics. 

3.  The  tip-path  plane  is  in  a  prescribed  position. 

4.  The  blade  dynamic  motions  can  be  described  by: 

a.  Three  symmetrical  bending  modes. 

b.  Teetering  motion,  and  two  antisymmetrical 
bending  modes. 

c.  Two  symmetrical  torsional  modes. 

d.  Two  antisymmetrical  torsion  modes  (with 

the  same  frequencies  as  the  symmetric  modes). 

5.  Control  is  applied  by: 

a.  One  symmetric  mode  and  one  antisymmetric 
control  mode,  each  giving  uniform  feathering 
along  the  blade  span. 
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b.  One  symmetric  mode  and  one  antisymmetric 

control  mode,  each  giving  differential  feathering 
of  the  inboard  and  outboard  sections  of  the  blade. 

When  the  second  and  all  higher  even  harmonic  root  shears  are  prescribed 
to  be  zero,  the  vertical  root  shear  transmitted  to  the  fu.c,elage  is  a  con¬ 
stant,  not  depending  on  azimuth  position.  The  computational  procedure 
can,  in  principle,  be  used  for  finding  the  pitch  inputs  required  to  elim¬ 
inate  canceling  root  shears  which  occur  at  odd  harmonics  for  the  teetering 
rotor,  but  numerical  results  for  this  case  are  not  presented. 

The  aerodynamic  forces  to  be  used  in  the  computational  procedure 
of  this  section  are  those  discussed  in  Section  1  without  the  simplifying 
assumptions  made  in  Section  2.  Thus,  reversed  flow  and  stall  effects 
are  to  be  taken  into  account  by  the  methods  described  in  Reference  6. 

DESCRIPTION  OF  ROTOR  DYNAMIC  MOTIONS 


The  dynamic  rn  i«~n  of  the  rotor  is  described  by  the  modes 
indicated  schematically  "  *gu  'e  11  which  are  applicable  for  the  case 
of  a  two-bladed  teetering  configuration.  All  of  these  modes  can. be 
handled  simultaneously  in  the  computer  program  which  has  been 
developed;  however,  in  general,  solutions  were  obtained  by  setting 
some  of  the  generalized  coordinates  for  these  modes  equal  to  zero. 

The  computational  procedure  which  was  developed  is  also  capable  of 
treating  the  case  of  the  three-bladed  articulated  rotor.  The  same  basic 
procedure  could  be  used  for  rotors  having  a  larger  number  of  blades  but 
would  require  some  modification  of  the  computer  program. 

It  can  be  seen  that  even  and  odd  subscripts  are  used  to  designate 
symmetric  and  antisymmetric  modes  respectively.  The  vertical  tip 
deflections  due  to  bending  in  the  first,  second,  and  third  symmetric 
bending  modes  are  denoted  by  ,  hv  ,  and  and  are  used  as  gener¬ 
alized  coordinates.  Similarly,  the  tip  deflections  h, ,  h}  ,  and  hF  are 
used  as  generalized  coordinates  for  the  teetering,  first  antisymmetric 
bending,  and  second  antisymmetric  bending  modes. 

It  will  be  noted  that  a  vertical  rotor  displacement  mode 
is  included  corresponding  to  vertical  motion  of  the  flapping  hinge 
produced  by  fuselage  and  hub  motion.  The  ho5  mode  is  introduced 
because  the  sum  of  the  generalized  aerodynamic  and  inertia  forces 
acting  in  this  degree  of  freedom  is  equal  in  magnitude  to  the  trans¬ 
mitted  shear.  However,  the  amplitude  of  the  hoS  motion  is  assumed 
to  be  zero  throughout  the  analysis. 


33 


VERTICAL  DEFLECTION  MODES 


SYH.  BLADE 
DISPLACEMENT 


hfti =  ANTISYM.  BLADE 
I  DISPLACEMENT 


(SAME  PITCH  CHANGE  (OPPOSITE  PITCH  CHANGE 

ON  BOTH  BLADES}  ON  OPPOSITE  BLADES} 


Co  -  1st  _YM. 
|  PITCH 


C|  =  1st  ANTISYM. 
|  PITCH 


Figure  11 


2nd  SYM. 
PITCH 


2nd  ANTISYM. 
PITCH 


SCHEMATIC  DIAGRAM  OF  MODES  USED  TO  DESCRIBE 
MOTIONS  OF  A  TWO-BLADED  TEETERING  ROTOR 
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An  antisymmetric  blade  displacement  >  ode  (hOA)  is  also 
shown  in  the  figure.  This  mode  has  no  physical  significance  in  the 
case  of  a  teetering  rotor,  and  motions  in  this  mode  are  assumed  to 
be  zero.  In  the  case  of  an  articulated  rotor  blade  with  an  offset  flapping 
hinge,  this  mode  corresponds  to  a  tilting  of  the  shaft.  The  generalized 
forces  acting  in  the  antisymmetric  h0A  mode  produce  canceling  root 
shears,  but  such  shears  can  transmit  moments  to  the  shaft  when  offset 
hinges  are  used. 

The  generalized  coordinates  @2 ,  @3,  and  6 y  used  for  the 
blade  torsion  modes  are  the  angular  deflections  about  the  elastic  axis 
at  the  blade  tip.  Since  the  pitching  of  the  blade  root  is  described  by 
the  control  modes,  the  angular  deflections  in  the  blade  torsion  modes 
are  assumed  to  be  zero  at  the  blade  root.  As  a  result,  the  deflection 
shapes  for  each  blade  are  the  same  for  the  symmetric  and  antisymmetric 
torsional  modes.  The  natural  frequencies  of  the  bending  and  torsional 
modes  indicated  on  Figure  11  are  presented  in  Appendix  I  for  the  con¬ 
figuration  studied. 

The  computational  model  permits  the  use  of  two  symme.ric 
and  two  antisymmetric  control  modes.  Ordinary  collective  pitch  is 
the  steady  component  of  the  motion  in  the  first  symmetric  control 
mode,  Ci  ,  while  ordinary  cyclic  pitch  is  the  first  harmonic  motion  in 
the  first  antisymmetric  control  mode,  c,  .  Higher  harmonic  motions 
can  be  introduced  in  both  the  c,  and  c2  modes  in  order  to  modify  the 
root  shears.  The  second  symmetric  and  antisymmetric  control  modes, 
and  Cj  ,  which  are  shown  schematically,  are  possible  differential 
control  motions  in  which  the  inner  and  outer  sections  of  the  blade  pitch 
in  opposite  directions.  The  mode  shapes  shown  for  c3  and  Cy  could  be 
replaced  by  other  radial  variations  within  the  framework  of  the  analysis. 
For  example,  a  linear  variation  in  pitch  might  be  used  which  could  be 
produced  by  a  device  introducing  a  moment  at  the  blade  tip. 

The  spanwise  variations  of  the  amplitudes  in  each  mode  are 
represented  by  dimensionless  f  functions.  Thus, 


■fy .  amplitude  of  motion  in  y*1*  mode  at 

'  ‘  radial  station,  r  ,  and  time,  t  . 


HARMONIC  AND  DISCRETE  TIME  TREATMENT  OF  VARIABLES 


The  investigation  has  been  limited  to  considering  steady  flight 
conditions  for  which  the  rotor  loads,  rotor  motions  and  transmitted 
shears  are  assumed  to  be  periodic.  However,  in  carrying  out  the 
aerodynamic  computations,  including  the  effects  of  the  wake  vorticity 
distribution,  it  is  convenient  to  utilize  the  values  of  the  variables  at 
discrete  time  increments  similar  to  the  treatment  used  in  Reference  6. 
As  discussed  in  Section  1,  the  positions  of  the  shed  and  trailing  vortices 
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in  the  wake  are  determined  from  uniform  inflow  theory  for.  the  given 
flight  condition.  Then,  when  the  values  of  Vt  ,  oe.  ,  a,  and  h  in  Equa¬ 
tion  (12)  are  given  at  the  A 'RA  blade  segment  positions  in  the  rotor 
aisk,  the  aerodynamic  problem  is  completely  defined.  The  solution 
to  this  aerodynamic  problem  gives  the  blade  circulation,  lift,  and 
pitching  moment  at  each  of  the  blade  segment  positions  which  is  equiva¬ 
lent  to  giving  the  variations  in  the  aerodynamic  loadings  at  discrete 
azimuth  or  time  increments.  The  transmitted  shears  and  generalized 
forces  acting  in  each  of  the  bending  and  torsional  modes  at  the  discrete 
azimuth  positions  can  be  computed  from  the  lift  and  moment  distributions. 
In  computing  the  dynamic  responses  to  periodic  forces,  however,  it  is 
desirable  to  use  z  harmonic  description  of  the  variables.  The  time 
variation  of  the  generalized  coordinate  can  be  approximated  by  a 
finite  Fourier  series  as  follows: 

M  =  X.  {(An  ii)  cos  n  #  +  (&n  ?i)  sin  n  t  ]  (3°) 

n  *  0 
f 

where  ~  (/\//j)  /I  ($i)t  > 

2.  (NA> 

(Ar>tc) 

z  (fJA) 

/  i  (M) 

cos(£-l)fr. 

The  symbol  (An$.c)  indicates  the  operation  of  multiplying  (fc)-A  ,  the  value 
of  the  generalized  coordinate  at  azimuth  position,  ^  =  Zrr(ib -1)/(HA),  by 
(2/(a/A)^co5  n^and  summing  for  (At A)  equally  spaced  azimuth  positions.  This 
is  simply  the  coefficient  of  the  nih  harmonic  cosine  term.  Similarly, 

(Bn  ^-)xs  the  coefficient  of  the  nih  harmonic  sine  term. 


DISCUSSION  OF  EQUATIONS  OF  MOTION  AND  ANALYTICAL  PROBLEM 

The  equation  of  motion  for  each  of  the  linear  deflection  and 
torsional  modes  can  be  written  in  the  general  form, 


-znG_ 


■(K.  Q.  -  CiT  Q. 


)l 


*V 


•> 


(32) 
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for  assumed  complex  periodic  variations  in  the  generalized  coordinates 
and  in  all  the  generalized  forces  acting  on  the  blade.  In  this  expression, 
the  are  generalized  mass  coefficients,  the  ^  ^  5  are  gyroscopic 

coupling^ coefficients,  the  5  are  generalized  stiffness  coefficients, 

the  are  centrifugal  force  coefficients,  and  the  <3‘}l  3  are  struc¬ 

tural  damping  coefficients.  Definitions  of  these  coefficients  are  given 
in  Appendix  I.  The  quantities  on  the  right-hand  side  arejfche  complex 
generalized  aerodynamic  force  acting  in  the  mode,  ;  the  gener¬ 

alized  centrifugal  force  arising  from  built-in  twist  and  preconing,  , 
the  complex  noncanceling  root  shear,  ;  and  the  complex  canceling 
root  shear,  Sc  .  The  shears,  S#  and  5C  ,  act  only  in  the  hoS  and  hon 
modes,  respectively,  as  indicated  by  the  <5  functions.  All  coefficients 
and  generalized  forces  are  the  values  per  blade. 

Complex  vectors  (indicated  by  bars)  have  been  used  in  Equa¬ 
tion  (32),  since  they  are  convenient  in  writing  the  conventional  expression 
for  structural  damping.  The  complex  periodic  variations  of  each  gener¬ 
alized  coordinate  (yt)  is  expressed  in  the  form, 

n  -0 


=  Z  [(A»fi)coS  %)sin  +L2-\ K Y-°^P  +  ( ^)sin 3 3 ) 


whex  ?  the  complex  coefficients  have  been  defined  in  such  a  manner  that 
the  real  part  of  Equation  (33)  is  identical  to  Equation  (30).  Similar 
expressions  hold  for  the  other  complex  periodic  quantities  ((?  ,Z  , 

>< 
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Only  the  real  part  of  Equation  (32)  is  considered  in  the  subse¬ 
quent  discussion  because  it  suffices  to  describe  the  real  pnysical  motion. 
After  substituting  expressions  of  the  type  shown  in  Equation  (33),  the 
real  part  of  Equation  (32)  can  be  arranged  in  the  form  of  a  sum  of  har¬ 
monic  terms,  and  the  coefficient  of  each  cosine  and  sine  harmonic  must 
vanish  for  the  complete  equation  to  be  satisfied.  The  requirements  on 
the  coefficients  of  cos  mp  and  sinnftin  the  equation  for  the  a**  mode  are 
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In  order  to  clarify  the  problem,  these  equations  and  certain 
auxiliary  constraint  equations  are  displayed  in  matrix  form  below. 


Form  of  equations  for  even  harmonics,  n  =  2,  4,  .  ,  . ,  12  (also  appli¬ 
cable  for  n  -  0  as  discussed  below): 
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It  follows  from  the  symmetry  of  the  two-bladed  teetering  rotor 
that  the  generalized  forces  acting  in  the  symmetric  modes  can  only  be 
of  even  or  zero  harmonic  order  while  those  acting  in  the  antisymmetric 
modes  can  only  be  of  odd  harmonic  order.  As  a  result,  the  response 
in  the  symmetric  modes  must  be  at  zero  or  even  harmonic  orders,  and 
the  response  in  the  antisymmetric  modes  must  be  at  odd  harmonics  as 
indicated  by  the  variables  appearing  on  the  left-hand  sides  of  Equa¬ 
tions  (36)  and  (37).  However,  this  does  not  mean  that  the  two  sets  of 
equations  are  decoupled,  because  coupling  can  take  place  between  har¬ 
monic  coefficients  of  all  orders  through  the  generalized  aerodynamic 
forces.  For  example,  the  combination  of  even  harmonic  blade  motion 
and  the  flow  dissymmetry  due  to  forward  flight  can  produce  odd  harmonic 
generalized  aerodynamic  forces. 

The  first  four  rows  in  Equation  (36)  represent  the  equilibrium 
of  the  harmonic  cosine  components  of  the  generalized  forces  in  the 
symmetric  vertical  deflection  modes  ( hoS  ,h ,  hb).  The  fifth  and 
sixth  rows  represent  the  equilibrium  of  the  cosine  components  of  the 
generalized  forces  in  the  symmetric  torsion  modes  ( Oz.Q^).  Similarly, 
the  ninth  to  the  fourteenth  rows  represent  the  equilibrium  of  the  sine 
components  of  the  generalized  forces  in  the  symmetric  vertical  deflection 
and  torsional  modes. 

Sixteen  variables  are  listed  in  the  column  matrix  on  the  left- 
hand  side  of  Equation  (36).  The  deflections  in  the  h0$  mode  are  assumed 
to  be  zero, but  the  noncanceling  harmonic  shear  coefficients,  (AnSfJ)  and 
which  are  associated  with  this  mode  are  treated  as  variables 
in  the  problem.  The  remaining  variables  are  the  sine  and  cosine  com¬ 
ponents  of  the  deflections  in  the  symmetric  bending  modes,  symmetric 
torsion  modes,  and  symmetric  control  modes.  There  are  no  equations 
of  motion  included  for  the  equilibrium  of  the  control  modes,  since  the 
characteristics  of  the  control  system  have  not  been  included  in  the  study. 
The  seventh,  eighth,  fifteenth,  and  sixteenth  rows  of  Equation  (36)  are 
four  constraint  equations  which  are  added  to  the  dynamic  equations  to 
make  the  number  of  equations  equal  to  the  number  of  unknowns  in  the 
column  matrix  on  the  left-hand  side.  Since  constraint  equations  are 
included  in  the  equation  set  for  each  harmonic,  a  well-defined  problem 
results  for  the  complete  system  of  equations  having  the  same  number  of 
equations  as  unknowns. 

The  implementation  of  a  pitch  control  system  for  the  elimination 
of  root  shears  might  be  based  on  several  different  concepts.  One  such 
system  might  apply  a  higher  harmonic  pitch  control  schedule  which  would 
be  selected  in  accordance  with  the  flight  condition  (i.  e.  ,  flight  velocity, 
rotor  rpm,  etc.  ).  Solutions  based  on  the  method  given  in  this  section 
indicate  the  pitch  control  schedules  which  should  be  used.  In  another 
possible  system,  the  output  of  a  sensor  measuring  the  blade  root  shears 
might  be  used  as  a  feedback  in  a  pitch  control  servo.  The  application 
of  such  a  servo  system  would  be  dependert  upon  obtaining  a  design  in 
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which  the  dynamic  motions  of  the  coupled  rotor-servo  system  would  be 
stable.  A  more  extensive  computer  program  would  be  required  to  study 
this  problem.including  the  dynamics  of  both  the  rotor  and  servo  systems. 

A  general  element  of  the  16  x  16  dynamic  response  matrix  on 
the  left-hand  side  of  Equation  (36)  is  denoted  by  D(l.J.n),  where  1  and 
J  indicate  the  row  and  column,  respectively,  and  n  the  applicable  har¬ 
monic  order.  The  [D]  matrix  has  been  partitioned,  and  the  forms  of 
typical  elements  in  the  submatrices  are  shown. 

The  harmonic  generalized  aerodynamic  forces  acting  in 
the  symmetric  modes  are  listed  in  the  first  column  matrix  on  the  right- 
hand  side  of  Equation  (36).  In  the  next  column  matrix,  a  list  is  given  of 
the  generalized  centrifugal  forces  (Z  In))  associated  with  preconing  and 
built-in  twist.  These  are  steady  forces  and  (n)is  zero  unless  n  =  0. 
The  last  column  on  the  right-hand  side  of  Equation  (36)  is  used  in  the 
computer  program  to  prescribe  the  values  of  four  variables  (for  n  -  0, 
only  two  of  the  cosine  coefficients  need  be  prescribed,  since  all  the  sine 
coefficients  are  zero). 

An  example  will  be  considered  in  order  to  clarify  the  meaning 
of  the  notation  used  in  the  constraint  equations.  If  it  were  desired  to 
use  the  seventh  equation  to  put  a  constraint  on  ( AnS (i.  e.  ,  the  n*h  har¬ 
monic  cosine  coefficient  of  noncanceling  root  shear),  a  setting, 

(38) 


would  be  used  in  the  inputs  to  the  computer  program.  This  means  the 
only  element  of  the  [£>]  matrix  in  the  seventh  row  (I  =  1)  would  be  in  the 
first  column  (J  -  l)  resulting  in  the  equation 

pp.  (An$n)  ~  pp  YP(n,7).  (39) 

The  input  YP(n,7)  would  be  set  equal  to  zero  if  a  solution  were  desired 
with  the  harmonic  noncanceling  root  shear  equal  to  zero. 

The  preceding  discussion  has  been  given  relative  to  the  even 
harmonic  equations,  but  analogous  considerations  hold  for  the  odd 
harmonic  equations.  Although  the  expressions  for  the  matrix  ele¬ 
ments  appear  to  be  the  same  in  Equations  (36)  and  (37),  different 
f  '*%■%/'>  Yf,  ?.■  >  %  pi'n<^Gyt7j  coe&icienls  are  involved, 

because  symmetric  and  antisymmetric  modes  occur  in  the  two  cases, 
respectively.  F^r  n-0  ,  there  are  no  sine  component  generalized 
forces  or  generalized  coordinates  and  somewhat  different  treatment 
is  required.  It  was  desired  to  use  matrices  of  the  same  order  for 
all  harmonics  for  convenience  in  the  computer  program.  The  n  =  0 
case  was  brought  into  this  framework  by  retaining  only  diagonal 
terms  in  the  ninth  to  sixteenth  rows  of  the  matrix  and  setting  the 
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right-hand  side  of  the  ninth  to  the  sixteenth  equations  equal  to  zero. 

Thus,  the  same  matrix  inversion  program  could  be  used  for  n-0  as 
for  the  other  harmonics. 

The  set  of  equations  of  motion  and  constraint  equations  for 
even  or  odd  harmonics  can  be  written  symbolically  in  the  form 

[D(l,Jtnj\[x(n,r)\  -  £  | GUI)]  +£*  \z(n,l)]  »  <40> 

where  appropriate  symbols  have  been  introduced  for  the  column  matrices 
in  Equations  (36)  and  (37).  The  [X(n,J)]  matrix  is  composed  of  variables 
corresponding  to  symmetric  modes  for  zero  and  even  harmonics  and  is 
composed  of  variables  corresponding  to  antisymmetric  modes  for  odd 
harmonics.  The  harmonic  generalized  aerodynamic  forces  are 
functions  of  the  lift  and  moments  on  the  blade,  which  depend  on  all  the 
harmonics  of  the  motion  as  mentioned  previously.  That  is, 

[&(n,  I)]  -  [f{(A0hz),  •  •  •  (A0C¥);  (A,h,)t  •  ■  •  (A,cs);  (3,  h,).--  (d,es  ); 

'  ‘  ’  >  ft  tz  hz  )j  '  ‘  '  (A/2  Cf) jj .  (41) 

An  explicit  expression  for  the  generalized  forces  in  the  forim 
of  Equation  (41)  is  not  available,  but  such  a  relationship  can  be  obtained 
from  the  solution  of  the  aerodynamic  problem  as  discussed  in  Section  1. 
For  given  harmonic  coefficients,  the  azimuthal  variations  of  the  gener¬ 
alized  coordinates  can  be  computed  using  Equation  (30).  The  input  vari¬ 
ables  or,«,and/i  in  Equation  (12)  can  then  be  expressed  in  terms  of  the 
generalized  coordinates  for  blade  motion  by  referring  to  Figure  4  and 
again  making  small  angle  approximations.  The  geometric  angle  of 
attack  at  each  blade  segment  position  L  relative  to  the  local  velocity 
la 

a*-lee*Z*cTCr}k  1421 

where  96  is  tne  built-in  twist  and  the  control  modes  are  as  shown  on 
Figure  11,  The  relat:ve  angular  velocity  due  to  the  motion  of  the  blade 
is 
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The  last  term  is  due  to  the  small  component  of  the  shaft,  angular  velocity 
along  the  blade  axis  and  was  neglected  in  most  of  the  computations.  The 
vertical  velocity  of  the  airfoil  relative  to  the  air  is 

At  -  jKt  s‘"  as  *  (vr  cos  )  4o<  <P)(a  *214 

■  *  L  K  \  ^ e?  %  l  ■  («> 

where  e  is  the  distance  between  the  midchord  and  elastic  axis  and  eQ 
is  the  distance  between  the  midchord  and  the  pitch  axis.  The  sum¬ 
mations  in  Equations  (42),  (43),  and  (44)  are  over  all  the  symmetric 
and  antisymmetric  modes  of  bending,  torsion,  or  control  motion. 

Having  determined  the  ak’s  and  hk ’s ,  Equation  (11)  can 
be  solved  for  the  bound  vortices  at  all  blade  segment  positions  in  the 
rotor  disk,  the  corresponding  Glauert  coefficients  can  be  evaluated, 
and  the  lift  and  pitching  moment  loadings  (/  and  17) )  can  be  computed 
by  Equations  (15)  and  (16).  It  is  then  a  simple  computation  to  deter¬ 
mine  the  generalized  forces  in  each  of  the  vertical  displacement  and 
torsional  modes  from  the  relations 

(&hi \  =  ft*  ( ^  M  dr  (45) 

and 

(Ge6)A  =  J -Pg  (r) dr ,  (46) 

where  the  subscript  £  indicates  the  generalized  forces  at  azimuth 
angle  The  harmonic  analysis  of  the  results  obtained  from  Equa¬ 
tions  (45)  and  (46)  is  the  final  step  in  relating  the  harmonic  coefficients 
of  the  generalized  aerodynamic  forces  and  the  harmonic  coefficients 
of  the  generalized  coordinates  of  the  blade  motions  as  indicated  in 
Equation  (41). 

Figure  12  is  a  schematic  diagram  of  the  analytical  problem 
which  must  be  solved  to  find  the  pitch  angle  inputs  necessary  to  elim¬ 
inate  oscillatory  root  shears.  The  purpose  of  this  diagram  is  to  show 
how  the  various  parts  of  the  problem  are  related  in  the  approach  which 
has  been  discussed.  The  overall  diagram  is  divided  into  two  boxes,  one 
showing  the  solution  for  the  aerodynamics  of  the  rotor -wake  system  and 
the  other,  the  solution  for  the  dynamic  blade  response,  required  pitch 
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Figure  12  SCHEMATIC  DIAGRAM  OF  THE  ANALYTICAL  PROBLEM  OF  FINDING 
PITCH  ANGLE  INPUTS  TO  ELIMINATE  OSCILLATORY  ROOT  SHEARS 
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angles,  and  residual  root  shears.  It  is  again  pointed  out  that  the  aero¬ 
dynamic  part  of  the  problem  is  treated  by  using  the  values  of  the  vari¬ 
ables  at  discrete  time  or  azimuth  angle  increments, while  the  solutions 
for  the  dynamics  motion  and  shears  are  treated  in  terms  of  their  har¬ 
monic  coefficients. 

ITERATIVE  SCHEME  OF  SOLUTION  FOR  HARMONIC  MOTIONS  AND 
- 

It  is  evident  from  the  schematic  diagram  shown  on  Figure  12 
that  the  dependencies  of  the  aerodynamic  loadings  on  the  dynamic  blade 
motions  are  very  involved.  The  loads  for  specified  input  motions  are 
obtained  only  after  an  iterative  solution  for  the  bound  vortices  as  dis¬ 
cussed  in  Section  1.  Consequently,  a  direct  solution  for  the  harmonic 
blade  motions  and  root  shears  due  to  the  aerodynamic  loadings  cannot 
be  carried  out,  and  an  iterative  scheme  of  solution  must  also  be  adopted 
for  their  determination. 

In  discussing  this  iterative  scheme  of  solution,  superscripts 
are  used  to  denote  the  values  of  the  variables  in  a  particular  iteration. 
Equation  (40)  can  then  be  written  as 


0(l.T,n) 


X%,J) 


n 


O(n.l) 


A 


Z(n,  r)J  + 


YP(nj) 


(47) 


where  |]  is  the  tth  approximation  for  the  column  matrix  of 

variables;  is  the  column  matrix  of  generalized  aerodynamic 

forces  based  on  the  preceding  approximation  for  the  column  matrix 
of  variables;  and  dl'\G-(r, ,  I)  ]  is  the  increment  in  the  column  matrix 
of  generalized  forces  from  the  (t-l)**  to  the  approximation. 

That  is, 


A(i)  \f(n,  I) J  =  \&  (t)(n,  Jj]  -  \g  (i  'f)(n,  I  j]  .  (48) 

No  superscripts  are  shown  for  ,  and  fYP]  matrices  because 

they  remain  constant  in  the  iterative  solution. 

A  possible  iterative  procedure  would  be  to  assime  AM[<s]zero 
in  computing  the  ^^approximation  for  the  column  matrix  of  variables, 
[Xft)]  ,  and  to  use  these  variables  to  compute  the  iih  approximation  for  the 
generalized  aerodynamic  forces,  [C-(t>]  ;  these,  in  turn,  are  used  to 
compute  [  neglecting  dttf,[G]  ,  etc.  However,  better  conver¬ 

gence  of  the  solution  would  be  expected  if  a  reasonable  estimate  could 
be  made  of  the  increments  in  the  generalized  aerodynamic  forces  from 
iteration  to  iteration.  A  relatively  simple  estimate  of  A(i)  [0]  can  be 
based  on  the  assumption  that  the  changes  in  the  aerodynamic  lifts  and 
pitching  moments  from  one  iteration  to  the  next  are  a  linear  function 
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of  the  corresponding  changes  in  the  harmonic  coefficients  of  the  blade 
motions.  Accordingly,  A(<'[iS]  is  approximated  by  an  expression  of 
the  following  form: 

A U) [G(n,I )  j  --  - [f  (I. I, «)][[* u)(nt /)] 

+  f  [H(l,T.n)]§X(t)(n+1,^-  \xli'f){n+1,  7)J} ,  (49) 

where  [f],  [f.]  ,  and  \H~\  are  square  matrices  which  are  functions  of  the 
harmonic  number,  n  .  Expressions  for  the  elements  of  these  matrices, 
which  are  based  on  quasi-static  aerodynamics,  are  presented  in  Ap¬ 
pendix  II. 

It  would  be  expected  that  the  increments  in  the  harmonic 
generalized  aerodynamic  forces,  would  be  primarily  caused 

by  the  changes  in  the  harmonic  blade  motions  and,  thus,  to  be  given 
by  the  first  term  on  the  right-hand  side  of  Equation  (49).  The  [Z]  and  [//] 
matrix  terms  give  the  coupling  of  with  variables  of  harmonic 

number  one  below  and  one  above  the  harmonic  number  of  the  generalized 
forces  under  consideration.  Couplings  of  this  type  would  be  expected  to 
be  of  first  order  in  the  advance  ratio,  //.  .  Couplings  with  harmonic  vari¬ 
ables  having  greater  separations  in  harmonic  number  (i,  e. ,  those  in  the 
matrices  (n-b2,T)\  ,  | X(n  -2,  J)] ,  etc.)  would  be  expected  to  be  of  second 
or  higher  order  in  the  advance  ratio  and  are  not  included  in  estimating 
A(t)[G(n,  Jj],  It  will  be  noted  that  the  [f]  and  [//]  matrices  produce  couplings 
between  the  generalized  r  rces  acting  in  symmetric  modes  and  the  motions 
in  antisymmetric  modes  and  vice  versa. 

After  substituting  Equation  (49)  into  Equation  (47),  the  following 
equation  results: 

[r(l,J,n)  ][xf<W)]  =  [YUl(n,J)]t  (50) 

where 

=  [Z?(X /,*)]  +  [£ (i,  J,  n)]  ,  (51) 
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A  m  \x(n  V- 1,  Jj]  =  | ~Xm(n  +  [x(i'%  V / ,  7)]  . 


(53) 


The  equations  obtained  by  setting  n  -  0,  1, ,  IE  in  Equa¬ 
tion  (50)  are  all  coupled  together  through  the[Z]  and  [#]  matrix  terms, 
giving  a  total  of  208  individual  simultaneous  equations  involving  the  208 
variables  which  are  the  elements  in  the  [X(n,J)\  matrices  with  harmonic 
numbers  ranging  from  0  to  12.  Although  these  are  linear  equations  when 
the  [ts]  matrices  are  considered  to  be  known  quantities,  it  is  imprac¬ 
ticable  to  obtain  a  solution  for  them  by  matrix  inversion  because  of  the 
large  number  of  variables  involved.  Consequently,  another  inner  loop 
iterative  solution  is  resorted  to  in  the  computational  procedure. 

The  computation  of  the  ^approximation  for  the  column  matrices 
of  harmonic  variables  starts  by  computing  the  [ Yq  ]  matrices  which  are 
the  [/'«]  matrices  of  Equation  (52)  neglecting  the  [L]  and  [//]  matrix  terms. 
This  decouples  the  equations  given  by  Equation  (50)  for  different  n's  and 
makes  possible  a  solution  of  each  matrix  equation  by  inverting  a  16  x  16 
[A]  matrix, 

[*"'■']-  HIV/"].  (54) 

The  prime  on  the  superscript  in  indicates  a  matrix  obtained  at  an 

intermediate  step  in  the  iterative  solution  of  Equation  (50).  Using  the 
^matrices  obtained  from  Equation  (54),  the  Ar<f![X]  matrices  in 
Equation  (53)  can  be  evaluated.  It  is  then  possible  to  use  Equation  (52) 
to  make  a  second  approximation  for  the  matrices  including  the  f/j 

and  \H)  matrix  terms.  This  inner  loop  iteration  process  can  be  repeated 
a  specified  number  of  times  (denoted  by  LGCD)  in  order  to  obtain  a  satis¬ 
factory  solution  of  Equation  (50). 

Once  the  [Y^]  matrices  have  been  obtained,  the  next  approxi¬ 
mation  for  the  generalized  aerodynamic  forces  can  be  computed.  Fig¬ 
ure  13  is  a  schematic  diagram  of  the  complete  computational  procedure 
showing  how  the  inner  loop  iterative  solutions  are  used  in  the  overall 
solution  for  the  harmonic  motions,  required  pitch  angles  and  residual 
root  shears.  The  overall  process  is  repeated  until  the  percentage 
variations  in  the  generalized  coordinates  from  one  overall  iteration  to 
the  next  are  below  prescribed  values. 
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HERAT  IOHS 


Figure  13  SCHEMATIC  DIAGRAM  OF  PRINCIPAL  STEPS  IN  ITERATIVE  SOLUTION 


4.  RESULTS  OF  COMPUTATIONS  FOR  PITCH  ANGLE  INPUTS  REQUIRED 
TO  ELIMINATE  NONCANCELING  OSCILLATORY  ROOT  SHEARS 


COMPUTATIONS  FOR  CONVENTIONAL  PITCH  CONTROL  (WITH  AND 

WITHOUT  "ELAT)  E-  T  OR'STON] - 


Flight  condition  DN67A  with  conventional  collective  and  cyclic 
pitch  controls  is  used  as  a  reference  case,  and  computations  for  this 
case  were  carried  out  to  compare  with  results  obtained  using  more 
complicated  control  inputs.  Two  sets  of  computations  were  made:  one 
including  deflections  in  the  torsional  modes  indicated  on  Figure  11 
(i.  e. ,  ,  &2  ,  63  ,  and  9<i)  and  the  other  setting  these  torsional  deflec¬ 

tions  equal  to  zero. 

The  cyclic  and  collective  pitch  angles  found  in  these  solutions 
are  presented  in  Table  III  and  are  the  ones  required  to  give  the  pre¬ 
scribed  thrust  per  blade  and  tip-path-plane  position  for  this  flight 
condition  (see  Table  I).  Plots  for  the  two  cases  were  superimposed, 
which  showed  the  azimuthal  variations  of  the  lift  loadings  and  bending 
moments  at  eight  radial  stations  along  the  span.  These  comparisons 
of  results  are  not  presented,  since  almost  identical  curves  were  obtained 
when  torsional  deflections  were  included  or  set  equal  to  zero. 

The  noncanceling  and  canceling  root  shears  obtained  in  these 
computer  runs  are  shown  in  Table  IV.  Here, again,  it  is  found  that  blade 
torsional  deflections  do  not  have  a  large  effect  upon  the  results  when  only 
zero  harmonic  collective  and  first  harmonic  cyclic  control  inputs  are 
used.  The  differences  in  the  amplitudes  of  the  harmonic  shears  computed 
with  and  without  blade  torsion  included  are  of  the  order  of  10  percent  of 
the  average  amplitudes  or  smaller. 

Convergence  difficulties  in  the  computations, which  were  associ¬ 
ated  with  torsional  deflections,  will  be  discussed  later.  These  difficulties 
were  overcome  in  computer  run  B-6, and  the  torsional  deflections  for  this 
case  were  found  to  be  small.  In  the  first  symmetric  torsional  mode,  they 
were  |0*|  =  0,813,  0.  164,  0.042,  and  0.  043  degrees  at  n  =  0,  2,  i,  and  6, 
respectively;  while  in  the  first  antisymmetric  torsional  mode,  they  were 
\&(\  =  0.  294,  0.  042,  0.  048,  and  0.  020  degrees  at  n  =  1,  3,  5,  and  7,  re¬ 
spectively.  Torsional  deflections  in  the  second  symmeti  .c  and  anti¬ 
symmetric  torsional  modes  were  an  order  of  magnitude  smaller. 

It  should  be  pointed  out  that,  in  the  calculations  made  including 
torsion,  the  blade  pitching  and  elastic  axes  were  assumed  coincident 
and  the  blade  was  assumed  mass  balanced  about  the  elastic  axis.  If 
the  blade  had  not  been  assumed  to  be  mass  balanced,  the  torsional 
degrees  of  freedom  would  probably  have  had  a  larger  influence  on  the 
results. 
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Table  Ml 

EFFECT  OF  BLADE  TORSION  ON  REQUIRED 
COLLECTIVE  AND  CYCLIC  PITCH  ANGLES 


FLIGHT  CONDITION  DN67A  WITH  CONVENTIONAL  CONTROL 


COMPUTER  RUN 

B-10 

B-6 

NO  BLADE 

INCLUDING 

TORSION 

TORSIONAL  MODES 

(A0CZ)  COLLECTIVE 

18.62  deg 

19.27  deg 

{A,  Ct)  CYCLIC  -  COS  COMP. 

3.56  deg 

3.26  deg 

(8,  Cf)  CYCLIC  -  SIN  COMP. 

-5.34  deg 

-5.36  deg 

CYCLIC,  TOTAL  AMPLITUDE 

8.42  deg 

6.27  deg 

CYCLIC,  PHASE  ANGLE 

56.30  deg 

58.70  deg 

Table  IV 

EFFECT  OF  BLADE  TORSION  ON  HARMONIC  ROOT  SHEARS 

FLIGHT  CONDITION  DN67A  WITH  CONVENTIONAL  CONTROL 


COMPUTER  RUN  8-10 


COMPUTER  RUN  B-6 


NO  BLADE  TORSION 


INCLUDING  TORSIONAL  MODES 


HARMONIC  ROOT  SHEARS  (ONE  BLADE)  * 


n 

AMPLITUDE 

lb 

PHASE  ANGLE 
(degrees) 

AMPLITUDE 

lb 

PHASE  ANGLE 
(degrees) 

0 

3200.0 

0 

3200.0 

0 

1 

2731.0 

-109.1 

2749.0 

-109.2 

2 

211.8 

170.4 

231.1 

169.7 

3 

253.3 

150.1 

236.5 

150.9 

45.5 

55.4 

45.2 

67.0 

89.5 

-43.4 

96.9 

-38.4 

5.0 

-141.9 

6.7 

-122.6 

7.6 

43.5 

6.5 

32.7 

8 

10.8 

128.3 

9.7 

135.6 

9 

2.0 

-64.9 

2.1 

-96.4 

10 

o79 

-3i.6 

10. 0 

-33.6 

mm 

2.5 

115.7 

2.6 

129.5 

m 

1.2 

-179.9 

l.l 

-180.0 

'ROOT  SHEARS  ARE  NONCANCELING  FOR  n  -  0  OR  n  =  even 
AND  ARE  CANCELING  FOR  n  =  odd 

NOTE:  AZIMUTHAL  VARIATION  OF  nth  HARMONIC  ROOT  SHEAR 
=  (Amp  1.)  x  cos  (iw/r+  Phase  Angle) 
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ELIMINATION  OF  NONCANCELING  ROOT  SHEARS  FOR  FLIGHT  CON- 
DllTON  DN67A  (BLADE  TORsIUETAl  DEFLECTIONS  AsSuMEff'zERQr 


In  this  section,  the  results  of  computations  are  discussed  which 
were  made  to  find  the  requirements  for  eliminating  noncanceling  root 
shears.  The  flexural  motions  of  the  blade  were  represented  by  the  first 
two  symmetric  and  first  two  antisymmetric  bending  modes, but  it  was 
assumed  that  there  were  no  deflections  in  the  torsional  modes.  Only 
the  first  symmetric  and  antisymmetric  pitch  modes  (cz  and  Ct)  were 
included, and  deflections  in  the  differential  pitch  control  modes  ( c 3  and  C+) 
were  set  equal  to  zero. 


As  mentioned  previously,  there  are  four  constraint  equations 
introduced  in  the  set  of  equations  for  each  harmonic  number  n  (see 
Equations  (36)  and  (37)),  and  four  of  the  harmonic  variables  for  each  n 
can  be  <?iven  prescribed  values.  The  procedure  used  for  making  the 
torsio..al  deflections  zero  was  equivalent  to  assuming  an  infinite  tor¬ 
sional  stiffness  and  did  not  involve  the  constraint  equations.  Two  of 
the  constraint  equations  at  each  n  were  used  to  make  all  the  harmonics 
of  C3  and  Cy  equal  to  zero. 


At  «  =  1,  the  remaining  two  prescribed  variables  are  the  cosine 
and  sine  components  of  the  teetering  motion  which  are  necessary  to 
establish  the  desired  tip-path-plane  position.  At  all  higher  odd  har¬ 
monics,  the  sine  and  cosine  components  of  the  inputs  in  the  first  anti¬ 
symmetric  control  modes  are  prescribed  to  be  zero; 


(A„cf)  --  (7  ( 

(&»  c ,)  -  0  J 


n  «  3,  S’,  7.9.  11. 


The  steady  component  of  the  root  shear  is  prescribed  to  be 
equal  to  the  desired  thrust  load  per  blade  of  3200  pounds.  The  con¬ 
straints  used  for  a  higher  harmonic  with  n  even  are  either  that  the 
harmonic  coefficients  of  the  root  shears  or  the  harmonic  coefficients 
of  the  first  control  mode  are  to  be  zero.  That  is, 


( -  0 

( )  =  O 
ov 

(  A„cz  )  -  o 


n  --  Z,  (,,8,  10,  iz. 


(dncz)  =  0  J 


Table  V  presents  data  on  the  harmonic  pitch  angle  inputs  re¬ 
quired  to  eliminate  noncanceling  root  shears.  The  first  two  columns 
repeat  the  data  given  previously  for  the  case  of  conventional  pitch  con¬ 
trol.  In  computer  run  B-12,  noncanceling  harmonic  shears  were  elim¬ 
inated  at  the  second  and  fourth  harmonics  by  applying  second  and  fourth 
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Table  V 

REQUIRED  PITCH  ANGLES  TO  ELIMINATE  NONCANCELING  ROOT  SHEARS 

FLIGHT  CONDITION  PN67A  (NO  TORSIONAL  DEFLECTIONS) 


■ 

COMPUTER 

RUN  B- 10 

COMPUTER 

RUH  B-12 

COMPUTER 

RUN  B-8 

mm 

PITCH  ANGLES  TO 

PITCH  ANGLES  TO 

BB 

PITCH  ANGLES  WITH 

ELIMINATE  ROOT 

ELIMINATE  ROOT  SHEARS 

CONVENTIONAL  CONTROL 

SHEARS  AT  n  =  2,4 

*T  n  =  2,4 

,6,8,10,12 

AMPLITUDE 

PHASE  ANGLE 

AMPLITUDE 

PHASE  ANGLE 

AMPLITUDE 

PHASE  ANGLE 

n 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

(degrees) 

(degrees ) 

0 

18.62 

0 

0 

I 

6.42 

56.3 

UllaS 

57.5 

2 

m 

1.21 

155.2 

.63 

41.7 

.62 

40.7 

.15 

-141.9 

1.01 

132.4 

10 

1  1 

1.10 

11.0 

12 

u  - - - 

.25 

-125.0 

AZIMUTHAL  VARIATION  OF  nth  HARMONIC  PITCH  ANGLE  =  (Antpl.)  x  cos  (n<7  +  Phase  Angle) 

Table  VI 

EFFECT  OF  HIGHER  HARMONIC  PITCH  ANGLE  INPUTS  ON  ROOT  SHEARS 


FLIGHT  CONDITION  DN6~*  (NO  TORSIONAL  DEFLECTIONS) 


COMPUTER  RUN  B- 10 

COMPUTER  RUH  8-12  j|  COMPUTER  RUN  B-8 

CONVENTIONAL 

CONTROL 

HIGHER  HARMONIC 
CONTROL  AT  n  =  2,4 

HIGHER  HARMONIC 
CONTROL  AT 
n  =  2,4,6,8,10,12 

HARMONIC  ROOT  SHEARS  (ONE  BLADE)  * 

AMPLITUDE 

PHASE  ANGLE  ! 

AMPLITUDE 

PHASE  ANGLE 

AMPLITUDE 

PHASE  ANGLE 

n 

(pounds) 

(degrees) 

(pounds) 

(degrees) 

(pounds) 

(degrees) 

0 

3200.0 

0 

3200.0 

0 

3200.0 

0 

1 

2731.0 

-109.1 

2756.0 

-109.1 

2756.0 

-109.1 

2 

211.8 

179.4 

o 

0 

0 

0 

3 

253.3 

150.1 

72.5 

144.6 

75.1 

154.6 

4 

45.5 

65.4 

0 

0 

0 

n 

5 

89.5 

-43.4 

66.5 

-48.2 

77.7 

-36.7 

6 

5.0 

-141.9 

5.3 

-104.4 

:  0 

0 

7 

7.6 

43.5 

8.5 

37,2 

18.3 

61.1 

8 

10.8 

128.3 

10.4 

131.8 

o 

0 

9 

2.0 

-64.9 

1.5 

-61.3 

17.4 

-66.0 

10 

8.9 

-31.6 

9.9 

-25.3 

0 

0 

11 

2.5 

115.7 

2.0 

101. 1 

11.4 

-  145.3 

!? 

_.K2- 

-179.9 

OA 

-1.9 

0 

0 

‘ROOT  SHEARS  ARE  NONCANCELING  FOR  n  =  0  OR  n  =  even  AND  ARE  CANCELING  FOR  n  =  odd 


AZIMUTHAL  VARIATION  OF  nth  HARMONIC  ROOT  SHEAR  =  (Amp I . )  x  cos  (n</'+  Phase  Angle) 
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harmonic  collective  pitch  control,  while  in  computer  run  B-8,  all 
noncanceling  harmonic  shears  were  eliminated  through  the  twelfth 
harmonic  by  applying  collective  pitch  control  at  all  even  harmonics 
through  the  twelfth.  The  root  shears  corresponding  to  these  pitch 
inputs  are  given  in  Table  VI  and  will  be  discussed  later.  Figure  14 
presents  a  comparison  of  the  pitch  control  schedules  corresponding 
to  the  harmonic  components  given  in  Table  V. 

The  pitch  angle  requirements  for  eliminating  the  second  and 
fourth  harmonic  root  shears  are  found  to  be  approximately  the  same 
whether  or  not  the  noncanceling  shears  for  the  sixth  through  the  twelfth 
harmonic  are  eliminated.  Noncanceling  shears  can  be  eliminated  at 
the  sixth  and  twelfth  harmonics  with  comparatively  small  pitch  inputs. 
However,  it  is  found  that  large  collective  pitch  inputs  are  required  at 
the  eighth  and  tenth  harmonics  if  the  root  shears  at  these  harmonics 
are  to  be  made  zero.  The  presence  of  these  large  higher  harmonic 
inputs  results  in  a  very  irregular  collective  pitch  schedule  as  shown 
by  the  curve  for  computer  run  B-8  on  Figure  14. 

A  physical  explanation,  of  the  need  for  the  large  higher  harmonic 
pitch  inputs  is  readily  obtained  by  considering  the  effect  of  the  dynamic 
response  of  the  rotor.  The  aerodynamic  forces  due  to  pitch  inputs  not 
only  produce  uot  shears  but  also  cause  teetering  and  bending  motions 
of  the  blade.  At  some  harmonics,  the  inertia  forces  resulting  from 
these  teetering  and  bending  motions  can  cause  root  shears  which  are 
nearly  equal  and  opposite  to  the  shears  due  to  the  aerodynamic  forces. 

As  a  result,  large  pitch  inputs  may  be  required  at  some  harmonics  to 
bring  about  small  changes  in  the  corresponding  root  shears. 

For  example,  when  conventional  control  is  assumed  at  flight 
condition  DN67A,  the  eighth  harmonic  root  shear  due  to  aerodynamic 
forces  has  an  amplitude  of  29.  2  pounds,  while  the  total  root  shear, 
including  dynamic  effects,  has  an  amplitude  of  10.8  pounds.  If  an 
eighth  harmonic  pitch  input  is  introduced  to  make  the  total  eighth  har¬ 
monic  root  shear  zero,  the  contribution  of  the  aerodynamic  forces  to 
the  root  shear  is  increased  to  an  amplitude  of  469.  2  pounds.  A  com¬ 
pensating  increase  in  the  inertia  forces  occurs,  making  the  resultant 
shear  zero. 

Thus,  for  this  case,  the  eighth  harmonic  transmitted  shear 
has  been  made  zero  only  at  the  expense  of  introducing  large  eighth 
harmonic  pitch  inputs  and  large  eighth  harmonic  aerodynamic  and 
inertia  force  loadings.  It  is  possible  that  this  situation  might  be 
improved  by  the  use  of  a  different  control  mode  such  as  differential 
rotation  of  the  inner  and  outer  sections  of  the  blade  which  might  change 
the  relative  magnitudes  of  the  aerodynamic  and  inertia  root  shears. 
However,  the  results  suggest  that  if  the  usual  control  mode  is  used 
which  gives  the  same  pitch  inputs  at  all  radial  stations,  it  might  be 
desirable  to  eliminate  noncanceling  root  shears  only  through  the  fourth 
harmonic. 
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f,  AZIMUTH  ANGLE  -  degrees 


Figure  14  COMPARISON  OF  PITCH  CONTROL  SCHEDULES  FOR  ELIMINATING 

NONCANCELING  HARMONIC  ROOT  SHEARS.  FLIGHT  CONDITION  DN67A 
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In  order  to  evaluate  the  penalties  resulting  from  limiting  the 
pitch  angle  inputs  in  this  manner,  the  corresponding  effects  on  the 
harmonics  of  the  transmitted  shears  are  shown  in  Table  VI.  It  is 
found  that  the  noncanceling  shears  at  n  =  6,  8,  10,  and  12,  which 
were  obtained  for  computer  run  B-12  (no  pitch  control  at  harmonics 
greater  than  n  -  4),  are  essentially  the  same  as  those  obtained  in  com¬ 
puter  run  B-10  (conventional  control).  The  canceling  root  shears  at 
the  seventh,  ninth  and  eleventh  harmonics  in  computer  run  B-12  are 
also  approximately  of  the  same  magnitude  as  those  obtained  with  con¬ 
ventional  control,  but  those  for  computer  run  B-8  are  considerably 
higher.  This  result  indicates  that  an  appreciable  interharmonic 
coupling  exists  through  the  aerodynamic  forces. 

The  azimuthal  variations  of  the  noncanceling  root  shears  corre¬ 
sponding  to  the  harmonic  data  in  Table  VI  are  plotted  on  Figure  15. 


5",  AZIMUTH  ANGLE  -  degrees 


Figure  15  AZIMUTHAL  VARIATION  OF  NONCANCELING  ROOT  SHEARS. 
FLIGHT  COHO  I Tl ON  DN67A 


Results  are  not  shown  for  computer  run  B-8  in  which  all  oscillatory 
noncanceling  root  shears  were  eliminated.  The  curve  for  this  run 
would  be  a  constant  value  of  -  3200  pounds  and  is  not  presented  to 
avoid  obscuring  the  other  data.  It  can  be  seen  that  when  pitch  control 
is  used  to  make  the  second  and  fourth  harmonics  of  the  root  shear 
zero,  the  principal  part  of  the  oscillatory  shear  is  eliminated. 
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Figures  16  and  17  show  the  effect  of  higher  harmonic  pitch 
control  on  the  oscillatory  lift  loadings  at  several  radial  stations.  The 
steady  components  of  lift  are  not  shown  on  these  plots.  On  Figure  16, 
the  results  for  computer  run  B-12  (pitch  control  to  eliminate  root  shears 
at  n  =  2,  4)  are  compared  with  those  for  computer  run  B-10  (conventional 
control).  On  Figure  17,  the  results  for  computer  run  B-8  (all  noncan¬ 
celing  root  shears  eliminated)  are  compared  with  the  results  for  the 
case  of  conventional  control. 

It  is  apparent  from  Figures  16  and  17  that  the  elimination  of 
harmonic  root  shears  by  pitch  angle  inputs  in  mode  c2  (i.  e.  ,  equal 
angle  changes  at  all  radial  stations)  does  not  eliminate  the  azimuthal 
variation  in  the  lift  loadings.  Large  high-frequency  variations  in  lift 
loadings,  particularly  near  the  blade  tip,  were  obtained  for  the  case  in 
which  all  noncanceling  harmonic  shears  were  eliminated  (see  Figure  17). 
When  only  the  second  and  fourth  harmonic  shears  are  eliminated,  the 
magnitudes  of  the  loadings  do  not,  in  general,  get  larger  than  those  with 
conventional  control  (see  Figure  16).  However,  the  changes  in  the  har¬ 
monic  content  due  to  higher  harmonic  control  result  in  differences  in  the 
azimuthal  variations  of  the  loadings  for  the  two  cases. 

Table  VII  presents  the  harmonics  and  maximum  peak-to-peak 
amplitudes  of  the  bending  moments  determined  for  the  same  conditions 
as  considered  in  Figures  16  and  17.  The  most  noticeable  changes  from 
the  case  of  conventional  control  are  large  reductions  in  the  second  and 
third  harmonic  bending  moments  and  a  reduction  in  the  peak-to-peak 
moments  at  the  inboard  stations.  The  bending  moments  from  the  eighth 
through  the  twelfth  harmonic  are,  in  general,  larger  for  the  case  when 
all  noncanceling  shears  are  eliminated  than  in  the  case  when  only  the 
second  and  fourth  harmonic  shears  are  eliminated, as  would  be  expected 
from  the  lift  loading  curves.  These  moments  would  probably  have  been 
larger  compared  to  the  lower  harmonic  moments  if  additional  bending 
modes  had  been  included  in  the  analysis  which  were  closer  to  resonance 
with  the  exciting  forces  at  the  eighth  and  tenth  harmonics. 

ELIMINATION  OF  NONCANCELING  ROOT  SHEARS.  COMPARISON  OF 
REQUIREMENTS  FCR.  FLIGHT  CONDITIONS  LN66ATH5N'67'A,"'ANI)  DN663 


The  results  considered  in  the  preceding  section  have  all  been 
obtained  for  the  same  basic  flight  condition,  DN67A.  Computations 
were  also  made  to  find  the  pitch  angle  inputs  required  to  eliminate  all 
noncanceling  root  shears  for  flight  conditions  DN66A  and  DN68A.  The 
required  pitch  angles  for  *he  three  flight  conditions  are  compared  in 
Table  VIII.  The  differences  in  the  required  second  and  fourth  harmonic 
pitch  angle  inputs  are  not  large  for  the  rpm  difference  of  flight  conditions 
DN67A  and  DN68A,  but  there  are  significant  differences  in  the  required, 
inputs  for  the  sixth  to  the  tenth  harmonic.  Large  differences  m  the  re¬ 
quired  pitch  angle  inputs  are  found  at  all  harmonics  for  the  two  conditions 
differing  only  in  flight  velocity. 
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LIFT  -  lb/ft  LIFT  -  lb/ft  LIFT  -  Ib/ft 


LIFT  -  Ib/ft 


f,  AZIMUTH  ANGLE  -  degrees 


SYMBOL 


COMPUTER 

RUN 


B-10 

B-8 


CONVENTIONAL  CONTROL 


PITCH  CONTROL  TO 
ELIMINATE  HARMONI 
ROOT  SHEARS  AT  n 
2,  6.  8,  10,  1 


FT  COMPONENTS 


ARE  NOT  SHOWN  ON  PLOTS 


Figure  17  EFFECT  OF  HIGHER  HARMONIC  PITCH  CONTROL  ON  OSCILLATORY  LIFT 
LOADINGS.  FLIGHT  CONDITION  DN67A.  COMPUTER  RUNS  B-8  AND  b-10 


K>  II  O 


N 


COMPUTER 

RUH 


PEAK 
TO  PEAK 


BENDING  MOM.  (in. -lb)  AT  r/R=0. 15 


HARMONICS  AT  WHICH  ROOT 
SHEARS  MADE  ZERO: 


raa 

n  =  2,4,6,6,10,12 

NONE 

ebb 

1773 

1831 

1006 

984 

3324 

3287 

2490 

2628 

627 

695 

469 

157 
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881 

1273 

517 

621 

615 

69 

83 

1619 

1913 

760 

649 

1 19 

116 

19 

16 

164 

140 

55 

59 

64 

271 

8.3 

8.3 

53 

27 

18 

19 

30 

203 

4.0 

3.7 

19 

42 

9.5 

7.3 

4 

!3 

0.4 

0.4 

10010 

10798 

7973 

6696 

BENDING  MOM.  (in. -lb)  AT  r/R=0.60 


HARMONICS  AT  WHICH  ROOT 
SHEARS  MADE  ZERO: 


NONE 


924 

444 

452 

413 

226 

593 

68 

78 

30 

24 

15 

7.6 

1.7 


IE 


=  2, 

3 

n  =  2,4,6,8,10,12 

NONE 

1 

=  2, 

BENDING  MOM.  (in. -lb)  AT  r/ R=0.36 


HARMONICS  AT  WHICH  R00T 
SHEARS  MADE  ZERO: 


=  2,7]^  =  2,4,6,8,10,12 
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2618 

160 

523 

69 

765 

16 


BENDING  MOM.  (in. -lb)  A'  r/R=0.80 


HARMONICS  AT  WHICH  ROOT 
SHEARS  MADE  ZERO: 


=  2,4  h  =  2,4,6,8,10,12 


6210 


Table  VIII 

REQUIRED  PITCH  ANGLES  TO  ELIMINATE  ALL 
NONCANCELING  ROOT  SHEARS 

COMPARISON  OF  RESULTS  FOR  FLIGHT  CONDITIONS  DN66A,  DN67A,  AND  DN68A 


FLIGHT 

CONDITION 

COMPUTER  RUN  D-2 

COMPUTER  RUN  B-8 

COMPUTER  RUN  C-2 

DN66A 

V  =  92. 1  KNOTS 

SI  =  313  RPM 

DN67A 

V  =  III.  4  KNOTS 

SI  =  313  RPM 

DN68A 

V  =  1  1 1 . 4  KNOTS 

SI  =  324  RPM 

n 

PITCH  ANGLE 

PITCH  ANGLE 

PITCH  ANGLE 

AMPLITUDE 

(degrees) 

PHASE  ANGLE 

(degrees) 

[jgjng 

PHASE  ANGLE 
(degrees) 

AMPLITUDE 
(degrees ) 

PHASE  ANGLE 
(degrees) 

0 

mm 

0 

mm 

0 

18.00 

0 

1 

59.1 

57.5 

57.3 

2 

0 

i 

161.3 

■H 

55.2 

1.06 

153.8 

0 

4 

.24 

8.4 

.62 

40.7 

.65 

39.8 

0 

6 

.64 

-116.2 

.15 

-141.9 

.30 

-140.4 

( 

8 

.66 

111.7 

1.01 

132.4 

1.17 

98.9 

10 

1.94 

4.4 

1.10 

11.0 

.25 

-5.0 

12 

.80 

-123.9 

.25 

-125.0 

.21 

-123.0 

AZIMUTHAL  VARIATION  OF  nth  HARMONIC  PITCH  ANGLE  =  (Ampl.)  x  cos  (n^+  Phase  Angle) 


59 


CONVERGENCE  DIFFICULTIES  IN  COMPUTER  SOLUTION 


In  the  first  attempt  to  solve  for  the  ordinary  cyclic  and  collective 
pitch  control  inputs  for  trimmed  flight  (including  the  effects  of  bending 
and  torsional  deflections),  the  solution  failed  to  converge.  It  was  found 
that  this  difficulty  resulted  from  the  natural  frequency  of  the  first  anti¬ 
symmetric  torsional  mode  being  very  close  to  the  frequency  of  the  fifth 
harmonic  exciting  forces.  Near  this  resonance,  the  inertia  and  spring 
forces  are  nearly  balanced  and  the  amplitude  of  the  torsional  response 
is  largely  determined  by  the  effective  aerodynamic  damping  forces. 

Under  these  conditions,  the  convergence  of  the  computational  solution 
is  quite  sensitive  to  the  accuracy  of  the  linearized  aerodynamics  used 
in  predicting  the  variation  in  the  generalized  torsional  force  from  one 
iteration  to  the  next.  By  analyzing  the  data  from  the  solution  that  did 
not  converge,  it  was  possible  to  make  a  better  estimate  of  the  aero¬ 
dynamic  coefficients  for  the  torsional  mode  than  was  obtained  from  the 
quasi-static  theory  of  Appendix  II.  When  these  revised  coefficients  were 
used  in  the  computation,  satisfactory  convergence  resulted. 

Convergence  difficulties  were  also  encountered  in  some  of  the 
computations  made  with  higher  harmonic  pitch  control.  The  failure  for 
convergence  in  these  cases  is  believed  to  be  due  to  the  sensitivity  of  the 
solution  to  the  estimated  aerodynamic  coefficients  near  resonant  con¬ 
ditions  arising  from  the  coupling  of  the  control  modes  with  the  torsional 
and  bending  modes.  Computations  were  not  carried  out  to  determine  if 
convergence  could  also  be  obtained  in  these  cases  by  improving  the  esti¬ 
mates  of  the  linearized  aerodynamic  coefficient  in  the  [£j  matrix. 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  following  specific  conclusions  are  reached  on  the  basis  of 
this  investigation: 

1.  The  computations  of  the  geometric  angle  of  attack 
required  to  eliminate  all  oscillatory  air  loads  show 
that  the  radial  variations  of  the  required  angles  are 
not  the  same  for  different  flight  conditions.  Although 
some  of  the  radial  variations  are  fairly  simple,  u 
number  of  the  radial  plots  of  required  angles  on 
Figure  9  showed  a  pronounced  "S"  shape.  Thus,  it 
is  concluded  th”t  it  would  be  a  difficult  undertaking 

to  design  a  pitch  control  system  for  the  elimination 
of  all  oscillatory  air  loads. 

2.  It  is  found  that  the  dynamic  response  of  a  rotor  to 
higher  harmonic  pitch  inputs  has  an  important  effect 
on  the  pitch  angles  required  to  eliminate  oscillatory 
root  shears. 

3.  All  noncanceling  vertical  root  shears  of  the  teetering 
two-bladed  rotor  from  the  second  through  the  twelfth 
harmonic  can  be  eliminated  with  higher  harmonic 
pitch  angle  control  introduced  at  the  blade  root  (as 

in  conventional  pitch  control  systems).  The  largest 
part  of  the  oscillatory  transmitted  root  shears  can  be 
removed  by  only  using  pitch  angle  inputs  at  the  second 
and  fourth  harmonics. 

4.  The  amplitudes  and  phase  angles  of  the  pitch  inputs 
required  to  eliminate  all  noncanceling  root  shears 
change  appreciably  with  flight  condition. 

5.  It  appears  undesirable  to  attempt  to  eliminate  all 
noncanceling  harmonic  root  shears  using  a  control 
system  with  pitch  inputs  introduced  at  the  blade  root. 
Higher  harmonic  pitch  inputs,  introduced  to  eliminate 
transmitted  root  shears,  excite  dynamic  blade  motions, 
and  the  inertia  forces  due  to  these  motions  also  pro¬ 
duce  root  shears.  It  is  found  at  some  harmonic  num¬ 
bers  (depending  upon  the  blade  parameters  and  the 
flight  condition)  that  the  root  shears  due  to  these 
inertia  forces  are  almost  equal  and  opposite  to  the 
root  shears  due  to  the  corresponding  aerodynamic 
forces.  This  means  that  large  pitch  angle  inputs 

and  large  blade  dynamic  motions  may  result  from 
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trying  to  eliminate  comparatively  small  root  shears. 

It  appears  undesirable  to  introduce  control  inputs  at 
harmonic  numbers  where  such  conditions  exist. 

The  results  of  this  study,  which  show  that  appropriate  pitch 
angle  inputs  can  produce  significant  reductions  in  the  transmitted 
oscillatory  root  shears,  offer  definite  encouragement  for  the  appli¬ 
cation  of  this  approach  to  reduce  fuselage  vibrations.  Howe-,  er,  a 
number  of  questions  remain  to  be  answer<  d  before  it  will  be  possible 
to  establish  the  value  and  feasibility  of  using  higher  harmonic  pitch 
control  for  this  purpose.  In  particular,  the  following  research  is 
recommended  for  answering  these  questions: 

1.  Determine  the  influence  of  the  pitch  angle  inputs 
which  are  introduced  to  remove  oscillatory  vertical 
forces  on  the  oscillatory  drag  and  side  forces  trans¬ 
mitted  to  the  fuselage.  {Some  effort  might  be  directed 
towards  adjusting  the  pitch  inputs  in  such  a  manner  as 
to  reduce  these  horizontal  forces  as  well  as  the  vertical 
forces.  However,  this  might  prove  to  be  a  more  difficult 
problem  because  of  the  nonlinear  variation  of  drag  forces 
with  angle  of  attack.  ) 

2.  Compute  the  blade  root  pi.ching  moments  on  a  rotor 
incorporating  higher  harmonic  pitch  control.  Deter¬ 
mine  if  the  reacting  moments  from  the  different 
blades  will  cancel  or  whether  they  will  tend  to 
excite  fuselage  vibrations. 

3.  Determine  how  helicopter  performance  is  affected 
by  the  higher  harmonic  pitch  inputs  which  are  intro¬ 
duced  to  eliminate  transmitted  shears. 

4.  Since  it  appears  that  there  are  large  differences  in 
the  higher  harmonic  pitch  angle  inputs  required  at 
different  flight  conditions,  it  is  recommended  that 

a  study  be  made  of  the  application  of  a  servo  system 
for  controlling  the  pitch  inputs  which  would  use 
measure'-'  transmitted  shears  as  a  feedback  signal. 

(Initial  study  of  the  stability  of  the  combined  rotor 
dynamics -servo  system  could  probably  best  be  done 
b-,  using  linearized  aerodynamics.  ) 

5.  All  computations  which  have  been  made  in  this  study 
have  been  for  the  case  of  control  inputs  applied  at  the 
blade  root  in  the  symmetric  or  collective  mode  (c?) 
and  in  the  antisymmetric  or  cyclic  mode  (c,)  .  How¬ 
ever,  other  modes  of  pitch  input  can  be  introduced 

in  the  computational  procedure  whicn  has  been 
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developed  (e.  g.  ,  symmetric  differential  motions 
of  the  inner  and  outer  sections  of  the  blade  (c¥)  or 
antisymmetric  differential  motions  (c3).  It  is  rpc 
ommended  that  the  application  of  pitch  modes  c3 
and  c *  be  studied,  particularly  at  those  harmonic 
numbers  where  inputs  in  the  usual  pi^ch  modes 
(ct  and  cz)  tend  to  excite  large  dynamic  responses 

6.  Finally,  it  is  recommended  that  a  study  be  made 
of  means  of  avoiding  the  convergence  difficulties 
in  the  computation  procedure  which  have  been 
encountered  at  resonant  conditions. 
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APPENDIX  I 

BLADE  PARAMETERS  USED  IN  COMPUTATIONS 


GEOMETRIC  PROPERTIES 


The  blade  planforms  consideied  in  the  computations  are  shown 
on  Figure  5.  A  preconing  angle  of  3.  0  degrees  is  used  in  all  compu¬ 
tations.  In  Section  4,  a  built-in  twist  (  &g)  is  included  which  varies 
linearly  from  zero  at  r  =  23  inches  to  -11.  0  degrees  at  the  blade  tip. 

NATURAL  FREQUENCIES  AND  BENDING  MODES 


The  natural  bending  frequencies  and  mode  shapes  for  the 
rotating  blade  are  based  on  computations  using  the  associated  matrix 
method  of  Reference  9,  while  the  torsional  frequencies  were  obtained 
using  a  Ralzer-type  computation.  A  summary  is  given  in  Table  IX  of 
the  frequencies  at  the  three  flight  conditions  for  which  computations  are 
carried  out. 


Table  IX 

NATURAL  BLADE  FREQUENCIES  USED  IN  COMPUTATIONS 


FLIGHT  CONDITIONS 
DN66A  AND  DN67A 

FLIGHT  CONDITION 
DN68A 

*>U 

(rad/sec) 

^9  C 

(rad/sec) 

FIRST  symmetrical  bending 

38,8 

1.18 

1 

1.18 

oJh  ,  SECOND  SYMMETRICAL  BENDING 

112.8 

3.43 

3.41 

t.Oh^  ,  THIRD  SYMMETRICAL  BENDING 

- 

- 

- 

- 

cJQ  ,  FIRST  SYMMETRICAL  TORSION 

168.8 

5.14 

169.2 

5.00 

Cc!a  ,  SECOND  SYMMETRICAL  TORSION 

440.9 

13.40 

441.5 

13.01 

,  FIRST  ANTI  SYMMETRICAL  BENDING 

96.4 

2.94 

99.2 

2.93 

Uu  .  SECOND  ANT  1 SYMMETR 1  CAL  BENDING 

171.6 

5.23 

175.6 

5.18 

C0Qf  ,  FIRST  ANTI  SYMMETRICAL  TORSION 

168.8 

5.14 

169.2 

5.00 

0Je  ,  SECOND  ANTI  SYMMETRICAL  TORSION 

3 

440.9 

13.40 

441.5 

13.01 
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Frequencies  are  given  for  the  first  and  second  symmetric  and 
first  and  second  antisymmetric  bending  and  torsional  modes.  The  fre¬ 
quency  and  mode  shape  of  the  third  symmetric  bending  mode  were  not 
computed.  The  deflections  in  this  mode  were  assumed  to  be  zero  in 
the  numerical  computations  for  the  pitch  angles  required  to  eliminate 
root  shears. 

COEFFICIENTS  IN  EQUATIONS  OF  MOTION 

The  equations  of  motion  as  presented  in  Equations  (34)  and  (35) 
of  Section  3  are  in  a  general  form  applicable  to  any  blade.  Expressions 
for  the  coefficients  in  these  equations  are  readily  obtained  for  specific 
configurations  as  done,  for  example,  in  Reference  8.  The  following 
assumptions  which  are  made  in  the  computations  of  Section  4  result  in 
a  number  of  the  coefficients  in  the  general  equations  becoming  equal 
to  zero. 

1.  The  relatively  small  gyroscopic  coupling  terms 
are  neglected. 

2.  The  elastic  and  pitch  axes  are  assumed  to  be 
coincident  and  to  intersect  with  the  shaft  axis. 

3.  The  blade  is  assumed  to  be  mass  balanced  about 
the  elastic  axis. 

4.  It  is  assumed  there  is  no  elastic  coupling  between 
the  bending  and  torsional  modes. 

For  this  case,  the  expressions  for  the  mass,  centrifugal 
force,  and  stiffness  coefficients  take  the  form 

z 

Iz  (root  fitting)  , 


Mi, 


mdr  +■ 


dr / r .o 


rmdr 


dr  *■ 


"  a  \K  T“’\  "hh 


> 
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where  77  is  the  mass  per  foot,  <Jlx/dzis  the  moment  of  inertia  per  foot 
about  the  elastic  axis  due  to  the  horizontal  distribution  of  mass,  and 
ctljdx is  the  moment  of  inertia  per  foot  about  the  elastic  axis  due  to  the 
vertical  distribution  of  mass.  The  vertical  distribution  of  mass  is 
neglected  in  treating  the  bending  modes  except  for  the  contribution  of 
the  root  fitting. 

The  numerical  values  of  the  nonzero  coefficients  which  were 
used  in  the  computations  are  summarized  in  Tables  X  and  XI. 


Table  X 

MASS,  CENTRIFUGAL  FORCE,  AND  STIFFNESS  COEFFICIENTS 
FOR  TEETERING  AND  BENDING  MODES 


"Mi  ’Mi  \  I,; 

(lb  sec^/ft)  (lb  sec^/ft)  (lb  sec^/ft)  (  lb/ft) 


TEETERING  (h,)  1.972  +2.847  -  1.972  0 

1st  SYMMETRICAL  BEND-  (h2)  1.553  2.246  -  2.107  57.5 

1st  ANT  I  SYMMETRICAL  8EN0.(/i3)  1.359  -0.891  -10.612  1143 


2nd  SYMMETRICAL  BEHD.  (h+)  K060  -0.780  -10.191  2450 

2nd  ANTI  SYMMETRICAL  BEND.  (hs)  1.327  0.470  -27.579  9236 


MASS,  CENTRIFUGAL  FORCE,  AND  STIFFNESS  COEFFICIENTS 
FOR  TORSION  AND  CONTROL  MODES 


Ma  *>. 
dc  6l 

(lb  ft  sec2) 

(lb  ft  sec2) 

K°C  Qi 
(lb  ft 
per  rad  ) 

1st  AXTISYMMETRICAL  TORSION  (0, ) 

0.2179 

-0.2294 

5961 

1st  SYMMETRICAL  TORSION  (fl2) 

0.2179 

-0.2294 

5961 

2nd  ANTI  SYMMETRICAL  TORSION  ( d3 ) 

0.2466 

-1.646 

45156 

2nd  SYMMETRICAL  TORSION  (0f) 

0.2466 

—  1  *  646 

46156 

M6,C,  = 

0.2965  lb 

ft  sec2  Td/C/  = 

-0.2965  lb 

= 

0.2965 

^2  ^Z 

-0.2965 

M*3Ci  = 

-0.1373 

V/ = 

0.1373 

MeH.cz  ~ 

-0.1373 

\CZ  = 

0.1373 

The  structural  damping  coefficient  is  assumed  to  be  ^  =  0. 
all  bending  and  torsional  modes. 


APPENDIX  II 


EXPRESSIONS  FOR  ELEMENTS  IN  [E],  [L],  AND  [H]  MATRICES 


In  estimating  the  increments  in  the  generalized  aerodynamic 
forces  from  one  iteration  to  the  next  of  the  computational  procedure, 
the  following  expressions  are  used  for  the  quasi-static  lift  and  pitching 
moment  loadings: 


■*  -  Zrrpb  f  Vth,  +  >5) 

yrj  (about  mtdchord)  -  -Vt h,  -  ~  b*cc  j  (56) 


where  V1  ,  cc ,  a:  and  h  are  given  by  Equations  (17),  (42),  (43),  and  (44), 
respectively.  The  incremental  generalized  forces  due  to  these  loadings 
are  computed  by  ?.  procedure  similar  to  that  used  in  Reference  10  and 
can  be  expressed  in  the  form 


% 


■rH 


(t)„ 


+  flA, 


Wi 


M. 

*  9i 


+  fLA 


9iV 


.«) 

*  h 


+  {(nSf  '<}sCn  +  (n*cf.  cos  +  >  (57> 

J  .  />  \ 

when  terms  proportional  to  /  are  dropped.  The  symbol  A  preceding 

the  symbol  for  a  given  quantity  is  again  used  to  denote  the  change  in 
that  quantity  from  the  (6- !)th  to  the  (i)th  approximation  of  the  iterative 
solution. 


Only  one-half  of  the  modes  of  the  two-bladed  teetering  rotor 
need  be  included  in  each  of  the  summations  in  Equation  (57)  because 
the  terms  involving  the  other  modes  are  zero  from  symmetry  consid¬ 
erations.  Terms  which  are  independent  of  yu  are  grouped  together  in 
the  first  sum.  The  only  modes  to  be  included  in  this  group  are 
those  having  the  same  symmetry  as  the  y.  mode.  Terms  proportional 
to  A  are  put  in  the  second  sum.  For  this  sum,  the  only  modes  in¬ 
cluded  are  those  having  the  opposite  symmetry  to  . 
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In  the  case  of  periodic  motion,  each  variable  in  Equation  (57) 
can  be  expressed  in  the  form  of  Equation  (30)  leading  to  the  following 
equations  for  the  harmonic  coefficients  of  Au)0  : 


-XfK*K4 

yl(H  ,yh-'  *  l^i  i  “"fc-’i) 

+(t  f;  y,  h }  ^*>4  is?) 

These  equations  are  combined  into  the  following  matrix  form  for  con¬ 
venience  in  the  computer  solution: 

4, \y}G(n,l) J  ;-[E(l,J,n)\  [ tf'xfn,  /)] 


+  J.  n)][itt]Mn- /,  7)J  +  j[h(i,J,  /?)]  [/ W/,  J)]  •  (60) 

In  this  expression,  \_X.(n.  jf\  is  the  column  matrix  which  lists  the  harmonic 
cosine  and  sine  coefficients  of  the  root  shear  and  the  generalized  coor¬ 
dinates  representing  the  blade  motions.  (The  forms  for  even  and  odd  n 
are  given  explicitly  in  Equations  (36)  and  (37).)  The  matrices  \X  In- ! ,  J)\ 
and  [X(n-bf,  correspond  to  the  (n-  i)th  and  harmonics. 


The  forms  of  the  16  x  16  [£],  [  2  ]  and  [H  ]  matrices  are  indi¬ 
cated  below  by  writing  the  expressions  for  typical  terms  in  the  sub¬ 
matrices  obtained  by  partitioning.  In  each  matrix,  the  first  and  ninth 
columns  are  zero  because  the  root  shears  do  not  appear  in  Equations  (58) 
and  (59).  Also,  the  seventh,  eighth,  fifteenth  and  sixteenth  rows  are  zero 
since  the  corresponding  equations  are  constraint  equations  rather  than 
an  expression  of  the  equilibrium  of  generalized  forces  in  the  conliol  modes. 
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Expressions  are  given  below  for  the  AtS  and  C  coefficients 
in  terms  of  the  blade  parameters.  The  A  coefficients  appear  only  in 
the  \E(I,J,  nj\  matrix  and  the  two  modes  coupled  by  each  of  the  terms 
must  have  the  same  symmetry.  On  the  other  hand,  the  3  and  C  co¬ 
efficients  appear  only  in  the  [Z.  ]  and  [  H  ]  matrices  and  couple  modes 
of  opposite  symmetry. 

Ah.h.  =  -rrp / bZ-Fh.-Fu.dr 

0  J  v  &  J 

Ah.e;  =  rr/>Jb:iefhAg,  dr 
t  j  i/  t  j 

Ah.c.  =  TrfjbP-e0fhfccLr 

4  J  J  fc  J 

=  rrpJb2ef6Jh.dr 
Ae.sj  -  -V/>fbZ(eZ+  ~)fe.  f0.dr 

A°i  bJ  =  fcj  dr 

bh.h-  **  "/Al  hr  -Fh.  fhdr 

A  J  l J  QC.  l'  J 

5  i  -  ZTrpJb  li£(i+e)^}  r/\  A dr 

\ij  =  z’r'’Jb{^(l+e‘>)*z\r%f‘jdr 
\  h  ‘  -Z7rfiJb  (i-i7e)  A  A  dr 

-  -fij e(.v&) tr  A  A  dr 
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